THE VISCOUS SURFACE-INTERNAL WAVE PROBLEM: 
NONLINEAR RAYLEIGH-TAYLOR INSTABILITY 

YANJIN WANG AND IAN TICE 

Abstract. We consider the free boundary problem for two layers of immiscible, viscous, incom- 
pressible fluid in a uniform gravitational field, lying above a rigid bottom in a three-dimensional 
horizontally periodic setting. The effect of surface tension is either taken into account at both free 
boundaries or neglected at both. We are concerned with the Rayleigh- Taylor instability, so we 
assume that the upper fluid is heavier than the lower fluid. When the surface tension at the free 
internal interface is below a critical value, which we identify, we establish that the problem under 
consideration is nonlinearly unstable. 



1. Introduction 

1.1. Formulation of the problem in Eulerian coordinates. In this paper we study the viscous 
surface-internal wave problem, which concerns the dynamics of two layers of distinct, immiscible, 
viscous, incompressible fluid lying above a flat rigid bottom and below an atmosphere of constant 
pressure. We assume that a uniform gravitational field points in the direction of the rigid bottom. 
This is a free boundary problem since both the upper surface, where the upper fluid meets the 
atmosphere, and the internal interface, where the upper and lower fluids meet, are free to evolve in 
time with the motion of the fluids. We assume that the upper fluid is heavier than the lower fluid, 
which makes the problem susceptible to the so-called Rayleigh- Taylor instability. We will consider 
the nonlinear Rayleigh- Taylor instability both with and without taking into account the effect of 
surface tension on the free surfaces. 

We will assume that the two fluids occupy the moving domain which we take to be three- 
dimensional and horizontally periodic. One fluid (-1-) fills the upper domain 

(1-1) ^+{t) = {y G t2 X M I rj^{t,yi,y2) < ys < 1 + r/+(t, yi, ys)}, 

and the other fluid (— ) fills the lower domain 

(1.2) n_{t) = {yeT^xR\-b<y3< r,^{t,yi,y2)}. 

Here we have written the periodic horizontal cross-section as T^ = (27rLiT) x (27rL2T), where 
T = R/Z is the usual 1-torus and Li,L2 > are fixed constants. We assume that 6 > is the 
constant depth of the rigid bottom but that the surface functions rj± are unknowns in the problem. 
The surface r+(t) = {2/3 = l + 7]^{t, 2/1,^2)} is the moving upper boundary of r_(t) = {y^ = 

rj-{t,yi,y2)} is the moving internal interface between the two fluids, and = {y^ = —b} is the 
fixed lower boundary of 

The fluids are described by their velocity and pressure functions, which are given for each t > 
by ^^±(*)') '■ ^^±(*) — ^ and p±{t,-) : ^±{t) — )■ R, respectively. The fluid motion is governed by 
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(1.3) 



the incompressible Navier-Stokes equations: 

p±dtv± + p±v± ■ Vv± + div S{p±,v±) = -gp±e3 in il±{t) 

divv± = in i^±{t) 

dtV± = V3,± - vi,±dy^r]± - V2,±dy2rj± on r±{t) 

S{p+,v+)n+ = peU^ - a+H+n+ on r+(t) 

v+=V-, S{p+,V-^)n^ = S{p-,V-)n^ + a^H^n^ on r_(t) 

V- = on S;,. 

Here the positive constants p± denote the densities of the respective fluids, g > is the acceleration 
of gravity, 63 = (0,0,1)^, and —gp±es is the gravitational force. Throughout the paper we will 
write IpI := /9+ — Since we are interested in the Rayleigh-Taylor instability, we assume that 
the upper fluid is heavier than the lower fluid, i.e., 

(1.4) p^>p_^ Ipj > 0. 

The quantity S{p±,v±) := p±I — p±^{v±) is known as the viscous stress tensor, where p± are 
the viscosities of the respective fluids and we have written / for the 3x3 identity matrix and 
0{v±)ij = djVi^± +diVj^± for twice the velocity deformation tensor. We have written div S{p±,v±) 
for the vector with z*^ component djS{p±,v±)ij; an easy computation shows that if divu-t = 0, 
then div S{p±,v±) = Vp± — p±Au±. The constant pe is the atmospheric pressure, and we take 
o"± > to be the constant coefficients of surface tension. We will assume that either a± = or 
a± > 0, i.e. we do not allow only one free surface to experience the effect of surface tension. In 
this paper, we let V=k denote the horizontal gradient, div* denote the horizontal divergence and A* 
denote the horizontal Laplace operator. Then the unit normal of r-i-(t) (pointing up), n±, is given 
by 

and H-i-, twice the mean curvature of the surface T±(t), is given by the formula 



(1.6) H± = div. 



^1 + |V,7?±|2 



The third equation in (II. 3p is called the kinematic boundary condition since it implies that the free 
surfaces are advected with the fluids. Notice that on r_(t), the continuity of velocity implies that 
= V-, which means that it is the common value of v± that advects the interface. For a more 
physical description of the equations (jl.3p and the boundary conditions in (jl.3p , we refer to [24J . 

Note that the constant 1 appearing above in the definition of ^+{t) and r-(.(t) is the equilibrium 
height of the upper fluid, i.e. the height of a solution with v± = 0, ry± = 0, etc. It is not a loss 
of generality for us to assume this value is unity. Indeed, if we were to replace the 1 in J^+(t) and 
r_|_(i) by an arbitrary equilibrium height L3 > 0, then a standard scaling argument would allow us 
to rescale the coordinates and unknowns in such a way that L3 > would be replaced by 1. Of 
course, in doing so we would have to multiply Li, L2, p±, cr±, g and b by positive constants. In 
our above formulation of the problem we assume that this procedure has already been done. 

To complete the statement of the problem, we must specify initial conditions. We suppose that 
the initial surfaces r-i-(O) are given by the graphs of the functions ^^±(0) = r/o,±, which yield the 
open sets ri±(0) on which we specify the initial data for the velocity, ^±(0) = vq^± : ^^±(0) — )• M^. 
We will assume that 1 + 770,+ > > —b on T^ and that r?o,±;^o,± satisfy certain compatibility 
conditions (see the discussion before Theorem 14. 4p required for the local well-posedness of the 
problem. We assume that the initial surface functions satisfy the "zero-average" condition 

(1-7) / Vo,± = 0. 

7X2 
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It is not a loss of generality to assume this. Indeed, if ()1.7p is not satisfied, then we could use 
the Galilean invariance of the equations and the flatness of the rigid lower boundary to restore the 
zero-average condition (see the introduction of [23j for details). For sufficiently regular solutions to 
the problem, the condition ()1.7p persists in time, i.e. 



(1.8) / ?7±(t) = for t > 0. 

Jt2 

Indeed, from the equations (jl.SP o. (ll.SP o. and ()1.3P g we obtain 



(1.11) 



(1.12) 



(1.9) — r/_ = / dft]^ = ti_ • n_ = / div?;_ = 0, 
at J J ^2 Jr^it) Jn-(t) 

and 

(1.10) / 1]+ = f dtr^^ = f v+ ■n+ = [ divv+ + [ v_ ■ ri- = 0. 

at Jt2 Ji2 JT+{t) Jn+(t) J^-{t) 

To simplify the equations we introduce the indicator function x ^'^^ denote 

V = v+Xn+{t) + v-Xn-{t), P = P+Xn+{t) + P-Xn^(t), 
P = P+Xn+it) + P-Xn^{t), P = P+Xn+{t) + P-Xn-(t), 

on J7(t) = ri+(t) U r2_(i). We similarly denote quantities on r{t) := T+{t) U r_(t), etc. We define 
the modified pressure through p = p + pgys — Pe — P+9- Then the modified equations are 

pdtv + pv ■ Vv + Vp = fiAv in il.{t) 

divv = in il.{t) 

dtf] = - vidy^r] - V2dy^r] on r(t) 

{P+I - P+^{v+))n+ = p+gri+n+ - a+H+n+ on r+(i) 

= on r_(t) 

{p+I — ^+ID'(u+))n_ — {p-I — p-B){v-))n- = Ipj gij^U- + U-H-U- on r_(t) 
(f,r/) |t=o= {vo,riQ)- 

In this paper we shall always unify the notations by means of (II. lip to suppress the subscript "it" 
unless clarification is needed. 

1.2. Reformulation via a flattening coordinate transformation. The movement of the free 
boundaries r-i-(t) and the subsequent change of the domains ^±{t) create numerous mathematical 
difficulties. To circumvent these, as usual, we will transform the free boundary problem under 
consideration to a problem with a fixed domain and fixed boundary. Since we are interested in 
the stability and instability of the equilibrium state, we will use the equilibrium domain. We will 
not use the usual Lagrangian coordinate transformation as in |181 [1] , but rather utilize a special 
flattening coordinate transformation that we introduced in [23]. 
To this end, we define the fixed domain 

(1.13) = U 0_ with f)+ := {0 < X3 < 1}, := {-h < X3 < 0}, 

for which we have written the coordinates as x G il. We shall write S+ := {^3 = 1} for the upper 
boundary, S_ := {^3 = 0} for the internal interface and Sf, := {x^ = —5} for the lower boundary. 
Throughout the paper we will write E = S+ U S_. We think of r]± as a function on E± according 
to r/+ : (T^ X {1}) x M+ ^ M and ??_ : (T^ x {0}) x M+ ^ M, respectively, where M+ = (0,oo). 
We will transform the free boundary problem in il(t) to one in the fixed domain by using the 
unknown free surface functions r}±. For this we define 

(1.14) 77+ := V+rjj^ = Poisson extension of 77+ into T^ x {X3 < 1} 
and 

(1.15) f)- := V-T]- = specialized Poisson extension of r]- into T^ x M, 
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where V± are defined by ()A.4p and ()A.9p . We now encounter the first key problem of how to define 
an appropriate coordinate transformation to flatten the two free surfaces together. Since we only 
need to transform the third component of the spatial coordinate and keep the other two fixed, 
we can flatten the domain by a linear combination of the three boundary functions, as introduced 
by Beale in [2j. However, this would result in the discontinuity of the Jacobian matrix of the 
coordinate transformation, which would then lead to severe technical difficulties in the proof of the 
local well-posedness of the problem. In [23] , we overcame this difficulty by flattening the coordinate 
domain via the following special coordinate transformation, writing h = 1 + X3/6: 



(1.16) 



B {xi,X2,x^ + a;|(r/+ - (1 + 1/^)??-) + 6r?~) = Q+{t,x) = (yi, 2/2, 2/3) G ^+{t), 
Vt^3 x^ {xi,X2,x^ + 6f/_) = e_(t,a;) = (yi, 2/2, 2/3) G ^~{t). 



Note that B(S_|_,i) = r+(t), 0(S_,t) = r_(t) and 0_(-,t) |s5= Id In order to write down 
the equations in the new coordinate system, we compute 




(1.17) ve = 1 and ^ := (Ve 





Here the components in the matrix are 
(1.18) 



A+ = xl {difi+ - (1 + l/6)ai??_) + difi^h, B+ = xl (52??+ - (1 + 1/6)92??-) + 927?-6, 
J+ = 1 + 2x3(7?+ - (1 + 1/6)??-) + xi(a3r/+ - (1 + \lh)d-if]J) + ^_/6 + a3f/_6. 



= 9i??_6, = 92?/- 6, 

j_ = i + f)^/b + a3??-6, K = j-^. 

Notice that J = det VG is the Jacobian of the coordinate transformation. We also denote 

W = dtQ^K 

N = {-diT], -^2??, 1) 

(1.19) 

T = ei + dir]e3 for i = 1, 2 



djQi for i,j = 1, 2, 3. 



It is straightforward to check that, because of how we have defined ??+, the matrix A and the other 
terms defined above are all continuous across the internal interface S_. This is crucial for the 
whole analysis in our proof of local well-posedness of the problem in [23]. Also, we may remark 
that one can replace in the coordinate transformation (jl.l6p by some smooth function h{x3) 
with h{l) = 1 and h{0) = h'{0) = • • • = 0, to let be more regular across the interface S_. 
Note that if ?? is sufficiently small (in an appropriate Sobolev space), then the mapping is a 
diffeomorphism. This allows us to transform the problem to one in the fixed spatial domain for 
each t > 0. We will transform the problem in different ways according to whether we take into 
account the effect of surface tension. 

In the case with surface tension, i.e. (T± > 0, ?? will be in a higher regularity class than u. This 
allows us to use the clever idea introduced by Beale in [2j to transform the velocity field in a manner 
that preserves the divergence-free condition. We define the pressure p on by the composition 
p{t,x) = po Q(t,x), but we define the velocity u on Q according to Vi o Q{t,x) = K6ijUj(t, x) 
(equivalently, Ui{t,x) = JAjiVj o 0(t,x)). The advantages of this transform are twofold. First, u 
has divergence zero in 0, if and only if v has the same property in Q{t). Second, the right-hand side 
of (ll.l2P n is replaced simply by U3. In the new coordinates, the system (I1.12P with surface tension 
becomes the equations for the new unknowns {u,p,r]) in the following perturbation form, denoting 
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the interfacial jump as If} = /. 



+ U'3=0 



-/- 



U'3=o: 



(1.20) 



pdtu — fiAu + Vp = f 
div u = 
dtV = U3 



in il, 
in 
on S 
on Sj 



M=0, lpI-fiB{u)je3 = {lpjg7]^+a-A,i]^)e3-g^ on S_ 



u_ = 

{u,r]) \t=o= (uo,%), 



on Sf, 



where the nonhnear terms / = {f^,f'^,f^) and g± = {g}^., g'^, g"^) are given by 
f = -pJAji[dt{K9jk)uk - Wd3{Kejk)uk + Kdk{K9ji)ukUi] 



(1.21) 



(1.22) 



(1.23) 



(1.24) 



+ nAjkdkiAjl)dlUi + lxJAniAjkAjldkdl{K0nm)Um 

+ pJAniAjkAjldl{K0nm)dkUm + p.J AniAjkAjldk{Kenm)dlUm 

+ p{AjkAji - 5ki)dkdiUi + ((5fej - JAjiAjk)dkP, i = 1,2, 3; 

C/Y = P-JAmkdk{K0jl)uiJ\fjl^ + p{Amk0jl - SikSjl)dkUlJ\fjT^ 

+ pJAjkdkiK0mi)uiMjT^ + n{AjkOmi - 6jk6ii)dkUiAfjT^ 

- pdiUidiT] - pdiU2d2'n + d^uiMidiTj 

- pdiUidiTj - pd2Uid2r] + pdkUsMkdit], i = 1, 2, 

5+ = 2^i{A3kKe3i - 53khi)dkUi + 2^iA3kdk{Ke-ii)ui 

+ 2p{Aikdk{Keji)ui + AikKejidkUi){NjNi - 5j35i3)W\-'^ 
+ 2p{A3kdk{Ke3i)ui + A3kKe3idkUi){\M\-^ - 1) 
-a+((l + |V*7?|2)-i/2-l)A,r? 

+ a+(l + |V*r?p)~3/2((5^^)252^ ^ 2air?a2r?aia2r? + (asr?)^^!!?); 

-5- = M JArakdk{Kejl)uiMjTi + {ArakOjl - SikSjl) {pOkUij AfjTi 

+ M JAjkdk{Ke^i)uiMjX, + (^jfc0 kSii)ydkUi}MjTZ 

- Ii4 diUidiT] - li4 diU2d2'n + Ip-d^uil Afidii] 

- 1/4 diUidiT] - [[/uj d2Uid2'n + I/i^fcUs]] A49jr/, i = 1, 2, 



-g^ = 2{A3kKesi - 5skS3i) IpdkUij + 2 [^1 AskOkiKe^ihi 

Aikdk{Keji)ui + Afei^^i/ [//afcU/l)(AAjM - 5j35a)\M\ 

1) 



-2 



(1.25) + 2(M A3kdk{Ke3i)ui + AskKOsi [/x5fcU;l)(|AA|-2 

+ (7_((l + |V,7?|2)~l/2-l)A,7? 

- f7_(l + I V*7?|2)-3/2 ((51^)2^2^ + 2dir]d2r]dld27] + {d2V?dl7]). 

Note that the coordinate transformation (|1.16p guarantees that [nj = on S_. 

In the case without surface tension, i.e., a± = rj will be in the same regularity class as 
f , so we can not transform the velocity field as above to preserve the divergence-free condition. 
Rather, we define the velocity u and the pressure p on directly by the composition with 0: 
u{t, x) = V o 0(t, x) and p = po Q{t, x). Hence in the new coordinates, the system (I1.12P without 
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surface tension becomes the following equations for the new unknowns (u, p, rj) in the perturbation 
form: 

pdtu — pAu + Vp = in Q 

div u = G'^ in $7 

dtV — U3 = G^ on S 

(1.26) { {p+I - ^+D(u+))e3 = p+gr]+e3 + G^ on S+ 

H=0, lpI-pB{u)}es = MgTl-es-G^_ on S_ 

U- = on 

^ {u,r]) \t=o= (no,r?o), 

where the nonlinear terms are given by 

(1.27) G^ = pWdsu - pu ■ Vau + K^a - A)ti - (V^ - V)p, 

(1.28) = (div - div^)n, 

(1.29) Gl = {p- p+gv+){e-i - M+) + ^i+(Btxe3 - B^+^x+AA+), 

(1.30) Gl = -{M - M 9V-)ie3 - M.) - M^ues - BavM)} , 

(1.31) = -uidiT] - U2d2'n. 

We may also refer to [lOj for more explicit expressions of G*. 



1.3. Main results. The Rayleigh- Taylor instability arises when steady states of two fluid layers 
with different densities are accelerated in the direction toward the denser fluid. The instability is 
well-known since the classical work of Rayleigh |17] and of Taylor |20] . For a general discussion of 
the physics related to this topic, we refer to [H] and the references therein. 

The Rayleigh- Taylor instability has attracted much attention in the mathematics community due 
both to its physical importance and to the mathematical challenges it offers. Presently, the linear 
Rayleigh- Taylor instability is well understood, but the nonlinear instability is much less so. The 
standard procedure (see, for instance, Chandrasekhar's book [1]) for analyzing the linear instability 
is to try to construct "normal mode" (horizontal Fourier modes of a fixed spatial frequency that grow 
exponentially in time) solutions to the linearization of problem ()1.12p . For instance, the inviscid 
(m± = 0) version of our present problem can be analyzed with normal modes as in [1]. It can be 
shown that without surface tension all the spatial frequencies are unstable. With surface tension, 
there is a critical frequency l^^lc = \/9 [pl / C- so that the spatial frequencies ^ with < |.^| < |^|c 
are unstable, while the spatial frequencies with |^| > |^|c are stable. This in particular implies 
that the inviscid version of problem (jl.l2p . defined in a horizontally infinite setting (i.e., horizontal 
cross-section in place of T^), is linearly unstable. However, in the horizontally periodic setting, 
if the domain is sufficiently small so that 



then there will be no frequency belonging to the interval < |^| < |^|c. This implies that the 
inviscid version of problem ()1.12p should be linearly stable. The same result for (viscous and 
inviscid) compressible fluids was shown by Guo and Tice [SJ E] . 

As a part of our previous paper [23], we proved that the viscous surface-internal wave problem 
p.l2p is nonlinearly stable for (t_ > cTc and that the solutions decay to the equilibrium state at an 
exponential rate. In this paper, we will continue our study of ()1.12p by showing that the problem 
is nonlinearly unstable for cj_ < cxc. 

Our main results can be stated as follows. The first theorem establishes the nonlinear instability 
for the case with surface tension. 
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Theorem 1.1. Assume that IpJ > 0, C7+ > 0, and < cr_ < ac, where Oc is defined by (jl.32p . 



be defined by (j4.17p and another (stronger) norm 



II 00 



be defined by KT8\i . 



Let the norm 

There exist constants 6*0 > and C > such that for any sufficiently small < l < Oq, there exist 



solutions 



to ([1:20]) so that 



77^(0) 



< Cl, but 



00 



> 



(1.33) 

Here the escape time > is 
(1.34) 

where A > 0, defined by ()2.32p . is the sharp linear growth rate obtained in Section\^ 



A i 



In [23| we also proved the nonlinear stability for the problem ()1.26p without surface tension 
((T± = 0), provided that \p\ < 0, i.e. the lighter fluid is on top. Our next theorem is the analog of 
Theorem 11.11 for the case without surface tension in the case [pj > 0. 

Theorem 1.2. Assume that [pj > and a± = 0. Let the norm |||-|||q be defined by (j4.17p and 



another (stronger) norm 

^0 > and C > such that for any sufficiently small < l < 6q, there exist solutions 
(fr26]l so that 



QQ defined by ()4.19p (with N > 3 an integer). There exist constants 

u^{t) 



to 



(1.35) 

Here the escape time T^ > is 
(1.36) 



u'{0) 
?7^(0) 



< Ci, but 



00 



u'(r') 
7?'(r') 



> 



1 1 ^0 

A* i 



where A* > 0, defined by ()2.33p . is close to the sharp linear growth rate A obtained in Section\^ 

Remark 1.3. The solution is actually a triple {u,p, rf), and the quantity |||-|||oo employed in Theorems 
and is actually a norm for {u,p,r]). See Remark \4-i\ for more details on our slight abuse 
of notion and notation. Also, although we have suppressed it in the notation, the quantity |||-|||oo 
appearing in Theorem depends on an integer N > 3 that measures the regularity of solutions 
(u € H^^ , etc). The choice N = 3 is the minimal regularity framework in which we can close 
our nonlinear instability estimates. However, by letting N be arbitrarily large we find that no 
high-regularity smallness criterion can prevent instability. 

Remark 1.4. Ln both Theorems and \1.2\ the instability occurs in the |||-|||q norm, which is just 
the L^ norms of u and rj. This means that although our instability analysis works in a framework 
with some degree of regularity, the onset of instability occurs at the lowest level of regularity. Note 
that an easy refinement of our analysis would allow us to replace |||-|||q with the norm of rj^. 
This highlights the fact that the instability occurs at the internal interface. We have chosen not to 
pursue this refinement in order to make |||-|||q consistent with a quantity that appears in our linear 
analysis. 

Remark 1.5. Theorems \l.l\ and \1.2[ together with our results in f23], establish sharp nonlinear sta- 
bility criteria for the equilibrium state in the viscous surface-internal wave problem. We summarize 
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these and the rates of decay to equilibrium in the table below. 





M<o 


M = o 


M >o 


a± = 


nonlinear ly stable 
almost exponential decay 


lorallv wpll-Dosprl 


nonliTiparlv nnstal^lp 


< cr+ 

< (T_ < CTc 


nonlinear ly stable 
exponential decay 


nonlinear ly stable 
exponential decay 


nonlinear ly unstable 


< cr+ 
0"c = CT- 


nonlinear ly stable 
exponential decay 


nonlinear ly stable 
exponential decay 


locally well-posed 


< cr+ 
CJc < (T_ 


nonlinear ly stable 
exponential decay 


nonlinear ly stable 
exponential decay 


nonlinearly stable 
exponential decay 



Note that our table identifies two critical regimes: o"_ = o"c, (T+ > 0, [pj > 0; and a± = 0, 
IpI =0. In these critical regimes we know from [231 that the problem is locally well-posed, but it is 
not clear to us what the stability of the system should be. 

In general, the passage from linear instability to nonlinear instability is delicate for partial 
differential equations due to, for instance, severe nonlinearities with possible unbounded derivatives. 
In 1995, Guo and Strauss [7] developed a bootstrap instability framework to treat such an issue. 
The main idea of the Guo-Strauss approach is to show that on the time scale of the instability, 
the stronger Sobolev norm of the solution is actually bounded by the growth of its weaker norm 
(see also Lemma 1.1 of [6j for the abstract framework of this method). For our Rayleigh- Taylor 
problem the term IpJ grj- is what leads to the instability; since it is of low order, we are led to 
use the Guo-Strauss bootstrap framework here. However, the bootstrap lemma in [6] can not be 
applied directly, so in the proof of Theorems 11.11 and 11.21 we use a slight variant. For the sake of 
completeness, we record an abstract version of this method in the following theorem. 



Theorem 1.6. Assume the following. 

(1) L is a linear operator on a Banach space X with norm 



and there exists another norm 



such that 



(1.37) 



|||e*^y|||o < Clc'^vW 

for some Cl and A > 0. 
2) N = N{y) is some nonlinear operator on X, and there exist another norm 
constant Cn such that 



II 00 



and a 



(1.38) 



||iV(y)|| 1112/1^0 



for all y £ X with |||y|||oo < oo. 
(3) There is a growing mode yo G X with |||yo||lo = 1 ^''^d lyolioo ^ ^''^^ 



(1.39) 



e^'yo = e^'yo 



for some A > with 4 < A < A. 



(1.40) 



(1.41) 



(4) There exists a small constant 6 such that for any solution y{t) to the equation 

d 



dt 



y + Ly = N{y) 



with |||y(t)|||oo — ^ '^^^ i ^ [0)^]> there exists Cs > so that the following energy estimate 
holds: 

Ibmlo < Cs Wymilo + \ f \\y{s)ll, ds + Cs f \\y{s)ll, ds + Cs f \\y{s)ll ds 

^ Jo Jo Jo 

for allt G [0,r]. 
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Let y'' be a solution to (jl.40p with initial data y^{0) = tyo- Then there exists a sufficiently small 
(fixed) number 6q > 0, which depends explicitly on Cl, Cn, Cs, X, A, yo, 6 but is independent of t, 
so that for any sufficiently small t > 0, 

(1-42) ll|y^(milo > y > 0. 

where the escape time T' > is defined by 

(1.43) r:=ilog^. 

A L 

Theorem 11.61 can be proved by a slight modifications of the proof of Lemma 1.1 of [6], which 
also follows along the lines of our proof of Theorems 11.11 and 11.21 and hence we omit its proof. 
In fact, one may read Theorem 11.61 as a description of our strategy for proving Theorems 11.11 and 
II. 2[ Comparing with Lemma 1.1 of [6], we see that (jl.37p and (jl.39p are two relaxations of the 
corresponding conditions in that lemma. This is precisely the case we encounter when we try to use 
the Guo-Strauss bootstrap framework to prove nonlinear Rayleigh- Taylor instability. Another point 
we want to mention here is that in Theorem 11.61 we implicitly assume that the linear growing mode 
yo can be used as initial data for the nonlinear problem ()1.40p . However, in our nonlinear problems 
()1.20p and (11.26P defined on a domain with boundary, to guarantee the local well-posedness the 
initial data must satisfy certain compatibility conditions that the growing modes to the linearized 
problem would not satisfy. Hence, we need to employ an argument from [12] to use the linear 
growing mode to construct initial data for the nonlinear problem in Section 14.31 

We will end this introduction by reviewing some previous mathematical results closely related 
to our current paper. The stability and instability analysis for the linearized problem without 
viscosity is well understood [4j. However, there is no general theory that guarantees the passage 
from linear instability to nonlinear instability for partial differential equations, so the question of 
nonlinear instability was not immediately resolved by the linear analysis. For the inviscid Rayleigh- 
Taylor problem without surface tension, Ebin [5] proved the nonlinear ill-posedness of the problem 
for incompressible fluids, Guo and Tice [9] proved an analogous result for compressible fluids, 
and Hwang and Guo [TT] proved the nonlinear instability of the incompressible problem with a 
continuous density distribution. For the viscous Rayleigh- Taylor problem with or without surface 
tension, Guo and Tice [8] proved the linear instability for compressible fluids, and with surface 
tension Priiss and Simonett [16] proved nonlinear instability for incompressible fluids in an LP- 
setting by using Henry's instability theorem. The proof of ^16j is quite abstract, hence in this paper 
for the case with surface tension we provide an alternative proof to show the nonlinear instability 
in the natural energy space. However, to our best knowledge, the result is completely new for the 
case without surface tension. In forthcoming papers, we expect to apply the framework here to 
show the nonlinear Rayleigh- Taylor instability for viscous incompressible fluids with magnetic field, 
the linear instability of which was obtained by Wang |22j . 

Notation. In this paper, for any given domain, we write for the usual L^-based Sobolev space 
of order A; > 0. We define the piecewise Sobolev spaces H^{^),k > by 

(1.44) H^{n) = {nxn± e H^{n±)] with norm ||n||^ := \\u\\%^f^^^ := + \\u-fH>={n.y 

When A: = 0, H^{VL) = H^{VL) = L^{n). We let the Sobolev spaces H'{T,±) for s G R be equivalent 
to H^{T'^), with norm || • we write 

(1-45) ||??||^ := lhllH=(s) •= ll^+ll//»(s+) + \\V-\\hs(t;_)- 

We do not distinguish the notation of norms on the domain or on the boundary. Basically, when 
we write and ||c^p||fc we always mean that the space is H^{Q), and when we write ||(9^r/||fc 
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we always mean that the space is When there is potential for confusion, typically when 

trace estimates are employed, we will clarify as needed. Let us define 

(1.46) oH^i^) := {u G H^{n) \ u \j^^= 0} and qHI{Q) := {u G qH^{VL) | divu = 0}, 

with the obvious restriction that the last space is for vector-valued functions only. It is easy to 
see that u G oH^{^) if and only if ti G H^{Q) with [nj = on S_ and u = on Sf,. We will not 
need negative index spaces on Q except for (o-ff^(S^))*. In our analysis, we will occasionally abuse 
notation by writing || • for the norm in {qH^{Q))* . Here it is not the case that {oH^{Q,))* = H^^ 
because of boundary conditions; we employ this abuse of notation in order to have indexed sums 
of norms include terms like || • ||4Ar_2j+i for j = 2N + 1. Sometimes we use || • \\lpx to denote the 
norm of the space U'{0,T; X). We will sometimes extend the above abuse of notation by writing 
II ■ IIlp//'-- for k = —1 in a sum of norms; when we do this we actually mean || • \\LP{oH^{n))'' ■ 

We shall also employ the notational convention for derivatives from |10j . When using space- 
time differential multi-indices, we will write N^^™ = {a = (qq, ai, . . . , a-m)} to emphasize that the 
0— index term is related to temporal derivatives. For just spatial derivatives we write N"^. For 
a G N^"*"™" we write = d^°d^^ ■ ■ ■ . We define the parabolic counting of such multi-indices 
by writing |a| = 2ao + ai + • • • + Um- We will also write V*/ for the horizontal gradient of /, i.e. 
V*/ = difei + 82/62, while V/ will denote the usual full gradient. 

For a given norm || • || and an integer A; > 0, we introduce the following notation for sums of 
spatial derivatives: 

(1.47) ||V^/f := \\dyf and ||VVf := E H^^/H'- 

\a\<k \a\<k 

The convention we adopt in this notation is that V* refers to only "horizontal" spatial derivatives, 
while V refers to full spatial derivatives. For space-time derivatives we add bars to our notation: 

(1.48) l|V^7f:= I|5"/f and llVVf := E H^^^H'- 

|Q|<fc \a\<k 

We allow for composition of derivatives in this counting scheme in a natural way; for example, 

||v.vS/|| = ||v^v,/|| = ||v^V||. 

We will employ the Einstein convention of summing over repeated indices. Throughout the paper 
C > will denote a generic constant that can depend on the parameters of the problem, integers 
used in energy definitions, and but does not depend on the data, etc. We refer to such constants 
as "universal." Such constants are allowed to change from line to line. When a constant depends 
on a quantity z we will write C = C{z) = Cz to indicate this. We will employ the notation a 
to mean that a < Cb for a universal constant C > 0. To indicate some constants in some places so 
that they can be referred to later, we will denote them in particular by Ci, C2, etc. 

2. Linear instability 
In this section, we consider the linearized equations of ()1.20p and (jl.26p in the unified form 



pdfU 


- fiAu + Vp = 


in 


divu 


= 


in Q, 


dtr] = 


U3 


on S 


(P+I 


- ^+B(n+))e3 = (p+OT+ - cr+A^77+)e3 


on 11+ 


N = 


0, Ipl - pB{u)} 63 = M 9V~ + a-A,r]-)e3 


on S_ 


u„ = 





on 



^ {u,7])\t=o = (no,??o), 

where C7± > 0. 
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2.1. Growing mode solution. We want to construct a growing mode solution to (j2.ip . so we 

assume an ansatz 

(2.2) u{t,x) = w{x)e^\ p{t,x) = n{x)e^\ r]{t,x') = C{x')e^\ 

for some A > (the same in both where x' = {xi,X2)- Substituting this ansatz into (12. ID . we 



find that C = wsls- By this fact we can eliminate C from (|2.ip and then arrive at the following 
time-invariant system for w and vf: 

Xpw — fiAw + Vvf = 



(2.3) 



div w = 

(yf+J - ^+D(u'+))e3 = X~^{gp+W3^+ - fT+A*W3,+)e3 
H = 0, I*/ - ^^I]){w)je3 = X-^{g {pj + a-A,W3)e3 



in Q 
in 
on 

on S_ 
on S/, 



Since the domain is horizontally flat, we are free to make the further structural assumption that 
the x' dependence of w, tt is given as a Fourier mode e*^ for ^ = (^i, ^2) £ L^^Z x L^^Z. Together 
with the growing mode ansatz (|2.2p . this constitutes a "normal mode" ansatz (see [3]). We define 
the new unknowns if, O^il^^-K : (—6, 1) ^ M according to 

-ie{x3)e'-'<, «;3(x) = ^(x3)e'^'^ 



(2.4) wi{x) = -iv9(x3)e"' W2{x) 

For each fixed nonzero spatial frequency ^, we deduce from the equations 
for ip,6,7p,-K and A, denoting ' = d/dx^: 



tt{x) = 7r(x3)e*''''^. 
23]) a system of ODEs 



(2.5) 



A2p(^ + ^A(|C|V 


- if") - A^ivr = 


in 




1) 


aV + mA(|cP^- 


- e") - A^vr = 


in 


i-b, 


1) 


AVV' + ^A(iepv 


- i)") + A^' = 


in 


i-b, 


1) 





in 


i-b, 


1) 


/^+A(6V'+ - 95'+) 


= ^+A(e2V+ - ^V) = 


at 


X3 = 


= 1 


—2p^Xip'_^_ + A7r+ 




at 


X3 = 


= 1 


iv^i = = m -- 

[mA(6V'-¥'')1 = 


= 


at 


X3 = 


= 


[^A(6^ - 0')1 = 


at 


X3 = 


= 


[-2/xAV'' + Avrl = 


--glpU-cT-\i?i^ 


at 


X3 = 


= 


= 6^ = = 





at 


X3 = 





We may eliminate vr by multiplying the first and second equations by .^i , ^2 respectively, summing, 
and then using the fourth equation to solve for 

(2.6) x\^\\ = -{x^pij' + xpm'^i^' - r')). 

Eliminating vr from the third equation in ()2.5p . we then obtain the following fourth-order ODE for 
^P■. 



(2.7) 



Given a solution ip to ()2.7p . we can recover the function vr via ()2.6p . To recover 93 and 6, we 
first note that the problem (12. 5|) is invariant under simultaneous rotations of {^p,0) and (i^i,^2)- 
Indeed, it is easily verified that if i? G 50(2) is any rotation operator, then R((p,6), R{£,i,(,2) is 
also a solution with the same ip and A. Then, given any ^ we choose a rotation operator so that 
-^5^ = (I'^l) 0)- It is easy to check that when ^2 = 0, we may simply take 9 = in (12. Sp . Then from 



-AViei'V'-^") 


= ^A(ieiv-2ieiV'+r") 


in 


(-&,!) 


p+x{\e^++i^'i) 


= 


at 


X3 = l 


p+x{i,% - 3|e|2^; 


) = AV+V'V + V+ + 


at 


X3 = l 


m = m = 




at 


X3 = ^ 


ipx{\e^+r)\-- 


= 


at 


X3 = ^ 


\px{r - 3iei V) 


l = [AWl + 9MieiV-'T-|eiV 


at 


X3 = ^ 


= = 




at 


X3 = -b. 



12 



YANJIN WANG AND IAN TICE 



a solution V'l^i, A, = {\(,\,0) to (j2.7p we can recover through the equation + = 0- 

Then we return to a full solution to (|2.5|) for ^ by setting {(p, 6) = R^'^{ip\^\^Q). 

The problem (j2.7p can be viewed as an eigenvalue problem with the eigenvalue A. However, 
due to the fact that A appears both quadratically and linearly, this problem does not possess a 
standard variational structure. In [8], Guo and Tice developed a robust method to overcome such 
a difficulty and constructed growing mode solutions for viscous compressible fluids. In [22], Wang 
also applied this method to construct growing mode solutions for viscous incompressible fluids 
with magnetic field in a slightly different setting. The method was also employed by Jang and 
Tice [12j for the compressible Navier-Stokes-Poisson system. The growing mode solution to (j2.ip 
can be constructed in the same way, though the actual procedure is somewhat simpler, so we will 
only sketch the procedure for later use, such as for deriving some estimates. For a more thorough 
presentation of the idea, we refer to [8] . 

In order to restore the ability to use variational methods we artificially remove the linear de- 
pendence on A in ()2.7p by defining the modified viscosities jl = s^, where s > is an arbitrary 
parameter. This results in a family (s > 0) of modified problems: 

" -A2p(|e|2^ - i;") = V - + r") m i-b, 1) 

sn+{\C\^iP+ + %b'l) = atx3 = l 

sf^+ii^'^ - 3iei W) = av+^v + <?p+iei v+ + v+ at xs = i 

(2.8) { [^^1 = ^=0 atx3 = 

[s/xdeiV + V'")! =0 atx3 = 

Isfiir' - 3|C| V)l = lAVV'l + g M 1^1 V - ^-lei V at X3 = 

ip- = ip'_ = at X3 = —b. 

Note that for any fixed s > and ^, (\2.8h is a standard eigenvalue problem for — A^, which has a 
natural variational structure that allows us to use variational methods to construct solutions. For 
this, we define the energies 



1 



(2.9) ^ 



-b 



+ ^iei'(^+iei' + 9P+)mi)\' + ^iei'(cT-ieP - g M)\m\\ 



(2.10) J{i^;m) = ^ l\m^\i;\^ + \ijf)dx3, 



b 

which are both well-defined on the space oH'^{{—b, 1)) := {v G H'^{{—b, 1)) | ipi—b) = ip'{—b) = 0}. 
We define the set 

(2.11) C = {V G oH^i-b, 1)) I J(^; lel) = 1}. 
Then we want to find the smallest — A^ by minimizing 

(2.12) _A2(|e|;s)=a(|e|;s) := mf E{i;;\C\,s). 

It is standard to show that a minimizer of ()2.12p exists and that the minimizer satisfies Euler- 
Lagrange equations equivalent to (|2.8p . A standard bootstrap argument also shows that the mini- 
mizer is in H^{{-b,Q)) (resp. i?'=((0, 1))) for all A; > when restricted to (-6,0) (resp. (0,1)), 
and hence is smooth when restricted to either interval. 

To construct the growing mode we then need to clarify the sign of the infimum (|2.12p . For 
0"- > (Tc, we always have that -E('0; \^\,s) > for any nonzero frequency ^ G L'^^'L x Lg^Z and 
any £C, which then means that q(|^|;s) > 0. This suggests that no growing mode solution to 
(12. ip can be constructed when (T_ > cTc, and then that the system should be linearly stable. In 
[23j we proved the nonlinear stability of the problem (I1.12P in the case (T_ > ac- However, when 
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< o"_ < (Tc and < |^| < |^|c := ■\/\p\ gjo- (it is interpreted that when (T_ = this means 
< 1^1 < cx)), then it is easy to see that if s > is smah enough we can always find ip G C such 
that E{'ip; |^|, s) < 0. Hence in this case there is ^s(|C|) such that E{ips(\(,\); |^|, s) = a(|^|; s) < 0. 

We want to show that there is a fixed point s so that A(|^|, s) = s, which will then allow us to 
construct a solution to the original problem (j2.7p . To this end, we study the behavior of a{s) as a 
function of s > 0. 

Lemma 2.1. For any fixed |^| G (0, |^|c), let a{s) = a(|^|;s) be defined by (I2.12p . Then the 
following hold. 



1) Q G C^^^((0, oo)) n C'^((0, oo)) and a{s) is strictly increasing in 



s. 



(2) There exist constants Co,Ci,C2 > depending on p±, fj,±, a±, b, g, |^| so that 

(2.13) a{s) < -Co + sCi, 
and 

(2.14) «(s)>-^|e|+5C2. 

P- 

Proof. We decompose the energy ()2.9p as 

(2.15) i?(V; lei, s) = sE^{^- lei) + i?o(V'; lei), 

where 



(2.16) ^i(^;|e|):=^ 



b 

(2.17) i?o(^; lei) := + 5P+)lV'(i)P + \\i\\cT.\i? - g M)\m?- 

The decomposition (j2.15p keeps the same form as in Proposition 3.6 of [8], hence (1) follows in the 
same way. 

Now we prove (2). First, for (|2.13p . note that for any fixed < |.^| < \^\c there exists V'l^i such 
that Co = — i?o('0|5|; 1^1) > 0. Then we have that i?(V'|g[; \i\-,s) < —Cq + sCi for some Ci > 0, and 
hence (I2.13P holds. Next, for (I2.14p . observe that for any ip G C, we have 



(2.18) 

> .\MA (|%_|jivi^.x3)' > 

Since the other terms in the energy E are nonnegative, we have 

(2.19) a(.) >-^M^ + . infill (V^). 

We denote by C2 this latter positive infimum; then (j2.14p follows and we conclude our lemma. □ 
For any fixed |^| G (0, |^|c), we then define the open set 

(2.20) 5 = a-^((-oo,0)) C (0,00). 

Note that S is non-empty and allows us to define A(s) = y^— a(s) for s £ S. We next show that 
there is a unique fix point s £ S so that s = A(|^|; s). 

Lemma 2.2. There exists a unique s £ S so that A(|^|; s) = y^— a(|^|; s) > and 

(2.21) s = X{\^\;s). 
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Proof. By Lemma |2. 11 there exists s=k > such that 

(2.22) 5 = a-^((-oo,0)) = (0,s*). 

We define A = y/—a on S and define the function <I> : (0, s=k) — )• (0, oo) by 

(2.23) ^s) = s/X{\^\;s), 

which is continuous and strictly increasing in s. Moreover, lim^-i^o ^(s) = and hm^^^^ ^(s) = +oo. 
Hence there is unique s S (0, s=k) so that ^>(s) = 1, which gives (j2.2ip . □ 

Notice that in Lemma [2T2l for any fixed |^| G (0, |.^|c) the fixed point s = s(|^|) E 5 is unique, 
we may write uniquely A(|^|) within (0, |i^|c), while the corresponding solution to (j2.7p is written 
by -0151- We have the following behavior of A(|^|). 

Proposition 2.3. The function A : (0, |^|c) — ?• (0, oo) satisfies the bound 

(2.24) Xm < 'jM, 

Moreover, 

(2.25) lim Xm) =0 for a±>0, and lim A(|C|) = f/ cj_ > 0. 

Proof We shah use the fact that ^(^-i^i; |^|, A(|C|)) = -A^d^l) < 0. First, we have 

\^?9M l0lfl(O)P 

(2.26) mmii^i^i) < ^^^^^'^1^1^ ^' , 

where Ei{tp^^^) := Ei{^^^^; |^|, A(|^|)). Since ip\^\{—b) = 0, we have 

(2.27) V|5|(0) = y ^V|5|^ix3 < ^(^y jV'ieiPd^aJ = ^/=p (y ^2M_|^p|4|l'd^3 
which implies that 

(2.28) \!iim^ .jmE,i^,^. 

This together with fTM implies (fZMl) . 

Next, we derive the limit behaviors of A(|^|). First, by (|2.14p . we have 

(2.29) < A^dei) < ^^lel, 

P- 

which implies that lim|g|^o ^^(1^1) = foi' — 0- ^^e other hand, by (|2.18p . we have 

(2.30) |V,|^l(0)|2<_i_. 
Hence, we obtain 

,2.31) A'(iei) < M!!Mt^|v,,^,(o)|= < ^afaM-'-"^!'), 

2 



which implies that lim|g|_^|^|^ A(|^|) = for cj_ > 0. □ 
By Proposition 12.31 we can then define 

(2.32) < A := sup < oo. 

o<l€l<l?|c 

For cr_ > 0, there is only a finite number of spatial frequencies ^ G {L^^^L) x [L'^^'L) satisfying 
ICI < \i\c-, so the the largest growth rate A must be achieved when < C7_ < Uc- For (T_ = it is 
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not clear whether A is achieved. However, we can achieve a growth rate that is arbitrarily close to 
A, and so in particular A* is achieved, where 

(2.33) < A/2 < A,, < A. 

We may now construct a growing mode solution to the linearized problem (j2.ip . 

Theorem 2.4. Let A be defined by (12.32^ and A* be defined by (12.33p . Then the following hold. 

(1) Let < (J- < (Tc- Then there is a growing mode solution to ()2.ip so that 

(2.34) Hmn^ = e^%m\H,, ||n(t)||^. = e^*||n(0)||^., for any k>0. 

(2) Let cr_ = 0. Then there is a growing mode solution to (|2.ip so that 

(2.35) \\rj{t)\\k = e'^**||ry(0)||fc, ||n(t)||fc = e^**||n(0)||fc, for any k>0. 

Proof Let |^| > be so that A(|^|) = A for cj_ > or A(|^|) = A=^ for cr_ = 0. Let tp^^i be 
the corresponding solution to (j2.7p with A(|i^|), which is constructed above by the minimization 
problem (j2.12p . We then define a solution to (j2.5p as described immediately after (j2.7p . which then 
allows us to define u, p, and rj according to (j2.2p . (j2.4p . and r] = X~^us\y:- Then we have that 
u G oH^{^) n if*^(0), p G H''{Q), and G H^{T,) for any A; > and {u,p,r]) solve the linearized 
problem ([2fT|) . Moreover, n,?? satisfy ([23^ or ([2?^ . □ 

2.2. Sharp growth rate. In this subsection, we will show that A defined by (j2.32p is the sharp 
growth rate of arbitrary solutions to the linearized problem (j2.ip . Since the spectrum of the linear 
operator is complicated, it is hard to obtain the largest growth rate of the solution operator in "L^" 
in the usual way. Instead, motivated by [SI [12]) we use careful energy estimates to show that e^* is 
the sharp growth rate in a slightly weaker sense (cf. ()1.37p ). 
First, we have the following energy identity. 

Lemma 2.5. Let {u,p,r]) solve (j2.ip . then 

+ \ I fi\ndtu\'^ = 0. 



Proof. We differentiate (12. ip ^ in time, multiply the resulting equation by dfU and then integrate by 
parts over Q. By using the other conditions in (12. ip . we obtain ()2.36p . □ 

The next result allows us to estimate the energy in terms of A. 
Lemma 2.6. Let u G qH^{Q) with crV*U3 G L'^ijl), then we have the inequality 

]- j a+lV^u^l^ + p+g\us\^ + ^ a-\V^us\^ -lp}g\u3\^ 

Proof. Let / be the horizontal Fourier transform of /. By the Parseval theorem, we have 
t; I cr+|V*n3p + p+5(|n3p + ^ / cr__|V*U3p - M 5|u3p 



2 

(2.38) 

" 47r2 
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Noticing that for ^ = 0, 
(2.39) iisiO) = 



U3 = I divn = 



and then 



(2.40) ^3(1) = f U3= [ 

we may reduce ()2.38p to be 



divn = 0, 



(2.41) 



= 4^ X Yl {^(^+1^1' + ^+5)1^3(1)1' + l{a-\C\' - M 5)1^(0)1'} 



Now for any fixed £, ^ 0, writing 

(2.42) ui{xs) = -iip{^,X3), U2{x3) = -i9{^,X3), usix^) = t/ji^xs) 

we have £,i(p + + 4^' = 0- The right hand side of (|2.41|) and the constraint ^np + + V'' = are 
obviously invariant under simultaneous rotations of ^ and {^p,0), so without loss of generality we 
may assume that ^ = (|^|,0) with |^| > and 6 = 0. Noting that for |^| G (0, \S,\c), by the definition 
of I^M), we have 

i(a+|C|2 + p+g)\us{l)\' + i(a_|eP - M 5)l%(0)|' 



(2.43) 



i(a+|eP + P+9mi)\' + ^(a_|^P - M <7)IV'(0)P 
^ (^(V^; lei, A) - ^ A/x(4|e| V? + llePV' + 



'"'MX3 



1 / .9 w , A 



>— -A^J{i;;\^\)-- / ;,(4|^|^|^'|^ + ||^|> + V'"|^)dx3 



lein 

= -Y j\pM' + \^\')dx3-^ J%i4\i;f + mi;-ipf)dxs. 

Here in the inequahty above we have used the following variational characterization for A, which 
follows directly from the definitions (|2.12|) and (|2.32p . 

(2.44) E{i;; ICI, A) > -A^J{i^; |e|), for any < |^| < and any V G oH\{-b, 1)). 

Note that for |^| > |^|c the left hand side of ()2.43p is nonnegative, so ()2.43p holds trivially, and so 
we deduce that it holds for all 7^ ^ G L^^Z x L^^TL. 

Plugging ()2.43p into (j2.4ip and translating the resulting inequality back to the original notation 
for fixed ^, by the Fubini and Parseval theorems, we find that 

\\ cr+|V*M3p +P+5k3p + ^ I cr-|V*U3p - M5IU3P 

2 A^ /"^ A /"^ 

(2.45) - ~4^ ^ ^ t/ ^I^(^3)I^^^3 + T / lA'^k^){^'i)\<^xz 

lit/ — o tJ — o 

A2 /■ , „ A 



Phi'- - / ii\^u\\ 

This is (12.37P and we conclude our lemma. □ 
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Now we can prove our main result of the subsection. 
Theorem 2.7. Let {u,p,r]) solve (|2.ip . Then we have the following estimates for t > 0; 

(2.46) \\u{t)\\l + Wdtuml + f \Hs)\\l ds < Ce'^'Wuml 

Jo 

and 

(2.47) Mt)\\l + \\dtvml/2 + r \\dsri{s)\\i/, ds < Ce^*(||r?(0)||2 + ||n(0)||i). 

Jo 

Proof. Integrating the result of Lemma 12.51 in time from to t, by Lemma 12.6^ we get 

I [ p\dtu\' + ]- f f ^,\o^tu\' 

^ Jn Jo Jn 

(2.48) =Ko-^(^j^ cj+|V*n3p + p+(7h3p + ^ a.\V^U3\^ - M g\u3\^^ 



/i|lL»u|^ 



with 
(2.49) 

Ko = U [ P\dtum^ + [ a+|V,tX3(0)P + P+g\u3m^ + / ^~|V.n3(0)|2 - Ipj gl^m' 
^ \Jn JE+ JE_ 

To compactly rewrite the previous inequality, we denote the weighted norms by 

(2.50) Ml := [ p\u\' and \\u\\l ■.= \ f ^i\^}u\'' 

Jn ^ Jn 

which are equivalent to || • ||o and || • ||f respectively. We denote the corresponding inner-products 
by {•,•)*, etc. Then (12.48P reads as 

(2.51) IWdtuml + 1^ \\dMs)\\lds <Ko + y ||n(t)||2 + ^\\u{t)\\l. 

Integrating in time and using Cauchy's inequality, we may bound 

rt 

|2 _ A IL./nM|2 



\u{s)\\l^ ds. 



A||ii(i)||^^ = A||n(0)||^,+A / 2{u{s),dsu{s))..ds 

(2.52) f ^ 

<A||n(0)||L+ / \\dsuis)\\lds + A^ [ 
Jo Jo 

On the other hand 

(2.53) Adt\\uit)\\l = 2A{u{t),dtuit)),<\\dtuml + ^^\Hml 
Hence, combining ()2.52p -( r2.53p with (I2.51|) . we derive the differential inequality 

(2.54) dtlWrnl + \\u{t)\\l <K, + 2A (^Mt)\\l + \\u{s)\\l ds^ 

for Ki = 2Kq/A + 2||n(0)||^^. An apphcation of Gronwah's inequality to (|2.54p yields 

(2.55) \\u{t)\\l + 1^ \\u{s)\\l < e'^^ml + ^(e^^* " 1). 
Now plugging ([235]) and (p32|) into ([23T]) . we find that 

(2.56) jWdtuml + \\uit)\\l < i^i + Ah(t)||^ +2A^* \Hs)\\lds < e'^\2A\\uml + K,) . 
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By the trace theorem and using the equations (j2.ip . we have 

(2.57) K,,K^ < \\um\l + \\dMO)fo < Mml 

and so (l236]l imphes I^AE\i . 

To prove (j2.47p . we use (j2.46p together with the kinematic boundary condition dtf] = and the 
trace theorem to estimate 



t rt 



ll^t^(*)IIV2 + / l|5s??(s)|li/2f^'5 = 1^3(^)11^1/2(2) + / ||n3(s)||^i/2(s)C^S 

(2.58) 



<ll^3(t)||?+ f\\u3is)\\lds<e'^'\\uml 
Jo 

and then by (f238D . 

(2.59) Ut)\\o < WmWo + r \\dsri{s)\\ods < e^*(||r?(0)||o + ||n(0)||2). 

Jo 

Hence, ([238]) and ([239]) imply (ITiT]) . □ 

3. Nonlinear energy estimates 

3.1. Energy estimates with surface tension. In this subsection, we wih derive the nonlinear 
energy estimates for the system ([1.20p . For this, we define the energy and dissipation, the definitions 
of which rely on the linear energy identity of the homogeneous form of ([1.20p : 

d f 1 f . .r, 1 f , ,0 ,_ ,9 1 



, , dt\2 
(3.1) 



+ 11 ^|ID)n|2 = 0. 
2 Jn 



According to this energy identity and the structure of the equations (11.20p . we define the instanta- 
neous energy (a higher regularity version of what we used in [23]) as 

2 1 3 

(3.2) £:=Y, \\d'tu\\l.2i + E \\9'tP\\l-2i-i + Ml + ^ ll^t^llt2m/2 

1=0 e=o e=i 

and the corresponding dissipation rate as 

2 1 3 

^ •= E 1 1 f^t"" 1 1 4-2^+1 + E l|f^tPll4-2£ + lkllll/2 + 5Z 1 1 '^t'^ 1 1 4-2^+5/2 

(3 3) 

+ Wdlvfo + II5'pIIh-i/2(e+) + III5'p1II^-i/2(s_)- 

We recall that in these definitions we have abused notation by writing || • := || • 

We will now derive our a priori estimates. Throughout the rest of this subsection we will assume 

that E(t) < (5^ for some sufficiently small 5 > and for all t in the interval in which the solution 

is defined. We will implicitly allow 5 to be made smaller in each result, but we will reiterate the 

smallness of 6 in our main result. 

To begin with, in the following lemma we present the estimates of the nonlinear terms / and g, 

defined by ([L2n) - (fT:25]l . in terms of £ and V. 

Lemma 3.1. Let f and g be defined by ([1.2ip - ([1.25p . then we have 

1 

(3.4) {\\9ff\\l-2i-2 + \\dfg\\Ui-3/2) < C£^ 

£=0 
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and 

2 

(3.5) {\\dlf\\l-2i-i + \\dlg\\l^2e-i/2) < C8V. 

Proof. The full expressions that define / and g are rather complicated, and a full analysis of each 
term would be tedious. As such, we will identify only the principal terms in the expressions and 
provide details for only the most delicate estimates. The other terms (lower order terms) may be 
handled through straightforward modifications of the arguments presented here. 

The principal terms appearing in / and g may be identified by first examining the regularity 
of the terms appearing in 8 and P, as defined by (I3.2p ~ (l3.3p . and by appealing to Lemma lA.ll to 
compare the regularity of f) to ry. We find that, roughly speaking, the regularity of u (and its time 
derivatives) is same as dtf], one order higher than p and at least one order lower than f). Based on 
this, we identify the principal terms in / and g defined by (ll.21j) - p.25p as 

(3.6) / ~ V^fyu + V^ryVn + Vr? jiV'^u 
and 

(3.7) g+ ~ V*7?+Vn+ + V^r/+-u+, g^ ~ V*7?_ [[/iVti]] + Vlr]-U. 
We will prove only the most involved estimate, namely, 

(3.8) \\dffe, + \\dh\ti/2<£-^- 

The other estimates may be derived through somewhat easier arguments. By Lemma lA.lHA.4l we 
have 

\\d?f\\-i := < \\d!v'fj\\o\\uh + mV^Udtuh + \\V'mo\\d!u\\i 

+ llSt^V^f/llollVnlli + \\dtV^fj\\o\\dtVu\\i + \\V^fi\\i\\d^Vu\\o 

+ l|5'Vf?||o||V2n||i + \\dtVf]\\i\\dtVMo + l|Vj?||3||M'v2^ll{oi?i)' 

(3.9) < \\d!vk/Mi + \\dtvk/2\\dtu\\i + ||r?||5/2l|5M|i 

+ \\dhh/2\W\\2 + \\dtvh/2\\dtuh + l|r/||5/2l|9M|i 

+ ll^t^'?lll/2ll'"l|2 + \\dtVh/2\\dtu\\2 + \\v\\7/2\\fJ'dt^^uk^Hiy 

<V£V. 

Here, in the last inequality, we have used the estimate 

(3.10) Wfid^vMiomr ^ II^Mli + Il5'v^x+||^_v2(2+) + \\d^ i/iVni < Vv. 

Indeed, by Holder's inequality and the trace theorem, we obtain that for any if G qH^{Q,) and any 
hj = 1,2,3, 

(3.11) 

{fidtdidju,(p)^ = - / fid^diu ■ djipdx + n+ / dfdiU+-Lp{ez-ej)- / 8"^ l^diu} ■ Lp{ez ■ ej) 

^ ll-^t^^lll Mil + ll'9t^Vu+||^-l/2(s^)||99||^l/2(£_^) + \\d^ If^Vuj ||//-1/2(S_)||^||//1/2(S_) 

^ (ll^t^^lli + l|5N'u+||^_i/2(s+) + \\dt IfiVuj Wh-i/^s.)) \Mi ■ 
Taking the supremum over such (p with \\(p\\i < 1, we get (I3.10p . 
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A similar application of these lemmas, together with trace estimates, implies that 

\\dh-\\-i/2 < l|9f'V,7?_||o|| If^Vuj 11,^2(2) + ||atV,r?_||o|| IfidtVuj h^^j,) 
+ l|V*r?_||2|| [/iS^Vn] ||h-i/2(e) + l|9'v2,7_||o||n||i2(s) 
+ \\dtVh-\\o\\dtuh2^j,) + \\Vh-\\o\\dMLHJ:) 
^ \\dtV-hhh + \\dtV~-\\i\\dtu\\2 + h-llall [M^Vn] || 



+ lia^Vlbhlli + ||5ir?_||2||5in||i + ||r/_||2||5M 



1 



and similarly, -^^2 ^ "^^^ Hence, p.Sp follows. □ 

3.1.1. Energy evolution. We first derive the energy evolution of the pure temporal derivatives. 
Lemma 3.2. Let {u,p,ri) solve (jl.20p . Then we have the estimate 

j2i\\d'Mml + \\dfrjml)+ fj2\\dlu{s)\\lds 

(3.13) ^=0 ■'^ ^=0 

t ft 



<C£{0) + C [ £{s)V{s)ds + C [ \\r]{s)\\lds. 
Jo Jo 



Proof. We multiply ()1.20p ^ by u and then integrate by parts respectively in ^1+ and by using 
the other conditions in (jl.20p . together with duality and trace theory, we then obtain that for any 
e > 0, 



l_d 

2di 



p\u\^+ a+\V,ri+\^ + p+g\T]+\^+ [ a.\V,r]A+\[ 
= f -u- g-u+ Mgr]-U3 



<Ce{\\f\ti + \\9\\\/2 + Ml)+e\\u\\l 



Now we integrate (j3.14p in time from to t and bound the resulting terms using the Korn inequality 
and the Poincare inequality on S_, as well as the nonlinear estimates (j3.5p : taking e sufficiently 
small in the resulting bound, we may absorb the ||u||f term onto the left to find that 

(3.15) \\u{t)\\l + Mt)\\\+ f\\u{s)\\lds<Cm + C f E{s)V{s)ds + C f Ms)\\lds. 

Jo Jo Jo 

Now we apply 5| with I = 1,2 to (jl.20p ^ and integrate and estimate as above to find that 

(3.16) = j dlf-dlu- I dlg-dlu+ I Mgdlr^^dlu, 

< C{\\dif\U + \\dl9\U,2 + \\dlri\U,2)\\dluh 

< CeiWdffWl, + \\dfg\ty, + WdfuWl) + eWdfuWl 

Here in the last inequality we have used the fact that 

(3.17) \\dfv\\-i/2 = I|5f-'n3||^-i/2(s) < CWdfuh. 
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Then we integrate ()3.16p in time from to t and argue as above to deduce that for i = 1,2, 

(3.18) \\d'Mml + \\dfvml+ f\\dfu{s)\\Us<C£{0) + C f £{s)V{s) ds + C f \\dtMs)\\Us- 

Jo Jo Jo 

Consequently, chaining ()3.15p and (|3.18p together, we get (|3.13p . □ 
Next, we derive the energy evolution of the mixed horizontal space-time derivatives. 
Lemma 3.3. Let {u,p,rj) solve (jl.20p . then for any e > 0, there exists constant > such that 

{Wu{t)\\l + Wr^{t)\\l) + f \\d''^i')\\lds 

(3.19) l"!-"^ 

<Ce£iO) + Ce [ £{s)V{s)ds + e f V{s)ds + Cef \\r]{s)\\lds. 
Jo Jo Jo 

Proof. Since the boundaries S-t and J^b are flat, we are free to take temporal and horizontal deriva- 
tives of the equations (|1.20p without disrupting the structure of the boundary conditions. Applying 
for a E N"*^"*"^ with \a\ < 4, ao < 2 to (jl.20P ]^. multiplying the resulting equations by d"u, and 
then integrating by parts over fi, we find that (again using trace estimates) 

1 d 
2di 



(3.20) 



Jo. Jy, / 2 jj^ 

Jo. Jt. Jt,- 

< c{\\dyu + ll5"5ll-i/2 + l|5"^ll-i/2)l|5"^x||i. 



Notice that for \a\ < 4, oq < 2, the kinematic boundary condition, the trace theorem, and the 
usual Sobolev interpolation allow us to estimate 

\\d''v\\-l/2 < \\v\\7/2 + \\dtV\\l/2 = M7/2 + ^311^3/2(2) 

< CsiWvfo + \Ml) + eiMl + hill) < c,M\l + IWWD + eT^- 

We then integrate ()3.20p in time from to t, plug ()3.2ip into the resulting inequality, and use 
Korn's inequality, Poincare's inequality, Cauchy's inequality, and the estimate p.Sp to find that 

ft 



Y{\\d^uml + \\d"r^{t)\\l)+ I Y \\d"Hs)\\lds 



|a|<4 " |a|<4 

(3.22) ao<2 ao<2 

<C£{0) + C f £{s)V{s)ds + e fv{s)ds + Ce !\ms)\\1 + \Hs)\\l) ds . 
Jo Jo Jo 

We then obtain the estimate (j3.19p from the above inequality and the estimate ()3.13p of Lemma 

E21 □ 

3.1.2. Full energy estimates. Notice that the estimates (|3.19p of Lemma [3.31 onlv contain the "hori- 
zontal" part energy estimates. In this subsection we will improve the estimates to be the full energy 
estimates. 

Theorem 3.4. Let {u,p, -q) solve (|1.2Up . If £{t) < 5"^ for sufficiently small S for all t G [0, T], then 
we have that 

(3.23) £{t)+ [ V{s)ds<C£{0) + C [ Ms)\\lds 

Jo Jo 

fort€[0,T]. 



(3.24) 
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Proof. As in [23], we will use the elliptic regularity theory for the stationary Stokes problems and 
the structure of the equations (jl.20p to improve the estimates. We divide our proof into several 
steps. 

Step 1 - £ estimates 

First, we may apply Of with £ = 0, 1 to the equations (|1.20p to find 

' -fiAdfu + Vdfp = -pdl+^u + d^J in 

div dl-u = in 

Idlp+I - /i+D(afu+)]le3 = p+gdf7]+e3 - a+A^dfr]+ + dfgl on S+ 

Idfuj = 0, Idfpl - fj,B{dfu)}e3 = Ipj gdfrj.ea + a.A^dfr]. - dfg^_ on S_ 

dfU- =0 on Sfe. 

Applying the two-phase Stokes regularity theory in Lemma IA.7I with r = 4 — 2£ > 2 to the problem 
([BTM]) and using (fXTUj) and ([23]), we obtain 

\\dfnml-2i + \\dfpmt2e-i 

(3.25) < II "^"^^114-2^-2 + \\'^tfi^)\\l-2e-2 + II ^114-2^-3/2+2 + ll^i ^114-2^-3/2 

^ W'-'t "'ll4-2f-2 II f^t'/ 1 1 4-2^-3/2+2- 

Chaining (j3.25p and (j3.19p . we obtain 

2 1 

(3.26) {\\9f<-2e + I|9'^ll4-2m) + E 11^*^^114-2^-1 <£' + 2, 

where we have written compactly Z to denote the right hand side of (|3.19p . 

Next, employing the kinematic boundary conditions and the trace theorem, we deduce from 
that 

3 3 
E ll^t ^(*)ll4-2£+3/2 — E 11^* "^^3(^)11^4-2^+3/2(2) 

(3.27) 

<i2\\dtMmi-2ii^i)<£'+z. 

Here, when £ = 3 we have used the fact ||c?f it3||/f-i/2(2) ^ ll'9|'"llo since divO^^n = in il. 

We complete the energy part of the desired estimate by summing (|3.26p - (|3.27p and then absorb- 
ing the E'^ term (taking b small enough) onto the left side of the resulting inequality. This yields 
the estimate 

(3.28) f(t) < Z. 

Step 2 ~ T> estimates for u and Vp 

Now we turn to the dissipation part estimates. Notice that we have not yet derived at estimate 
of rj in terms of the dissipation, so we can not apply the two-phase elliptic estimates of Lemma [A.7l 
as above. It is crucial to observe that from ()3.19p we can get higher regularity estimates of u on 
the boundaries S = U Indeed, since S+ are flat, by the definition of Sobolev norm on 
and the trace theorem, we may deduce from ()3.19p that for ^ = 0, 1, 2, 



1^ l|5N(^)ll^4-2.+v2(s) ds < 1^ {\\dfu{s)\\l,,.^^^ + \\vf^dfu{s)\\l,,,^^^ 

j[*(l|5M-)ll? + l|Vr^^)5fn(.)||?) ds 



(3.29) 

< / fliaM5)ll^ + l|Vr^^^5fn(5)|in ds<^. 



ds 



This motivates us to use the one-phase Stokes regularity theory of Lemma IA.81 
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The key point is to observe that (dfu-^., dfp+) with £ = 0,1 solve the problem 

r -/i+Aa|n+ + vdfp+ = -p+dl^\+ + dff+ in n+ 

(3.30) < divdfu+ = in n+ 

[ dfu+ = dfu^ on S, 

and {dfu-,dfp^) solve the problem 



(3.31) 



div 9j u_ = in 
dfu- = dfu^ on 
dfu- = on Sfo. 



We apply Lemma fA.SI with r = 4 — 2£+ 1 to the problems (j3.30p and (j3.3ip respectively, using the 
estimates (j3.5p . (|3.29p . and then summing up, to have 



r \\dfu{s)\\l_,,^, + \\Vdfp(s)\\l_,,_,ds 
Jo 



(3.32) 



< \\dtMs)\\L2e-i + \\dfm\\Ue-i + l|9fn(.)||^,_,,+,/,(^) ds 

< f\\d',^Ms)\\L2e-ids+ f £{s)V{s)ds + Z 
Jo Jo 



< f\\d'+Ms)\\l-2e-ids + Z. 
Jo 

Chaining (|3.32p and (|3.19p . we obtain 

l-t 2 1 
(3.33) / 115^(^)111-2^+1 + E l|V5fp(«)llt2^-i ds < Z. 

e=o e=o 



Step 3 -T) estimates for r/ and p 

Now that we have obtained (j3.33p . we estimate the remaining parts in T). We will turn to the 
boundary conditions in ()1.20p . First we derive estimates for rj. For the time derivatives of t/, we 
use the kinematic boundary conditions, the trace theorem and (j3.33p to obtain 



/ Y\\9fr]{s)\\l_2i_^^/2ds = Y\\dt^U3{s)\\H4-2i+5/2(^s)ds 
Jo Jo 



(3.34) 

ft 'J 



^ I T.\\^t-'Ms)\\Uie-i)+ids<Z. 



=1 

For the term rj without temporal derivatives we use the dynamic boundary conditions 

(3.35) - cr+A*?7+ =p+- 2/^+53^3,+ - g+- p+gr]+ on S+ 
and 

(3.36) - fT_A*7/_ = l-p + 2^53113]] - g- + Ipj gri_ on S_. 

Notice that at this point we do not have any bound on p on the boundary S yet, but we have the 
estimate of Vp in Cl. Applying V^, to (j3.35p - (j3.36p and employing the standard elliptic theory on 
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as well as the trace theorem and (j3.33p . we find that 

\NMs)\\l/2ds < ||V*p(s)||^5/2(s) + II V.a3tX3(s) 11^5/2(2) + \\^*V{s)\\l,^ds 

(3-37) < j'^ \\V.v{s)\\l + ||«(s)||i + ||ry(.)||?/2 ds 

<Z + j\\ri{s)\\l^ds. 

Using Poincare's inequality on Ti± (since Jrp2 = 0) and the same type of Sobolev interpolation 
trick we used in (|3.2ip . we deduce from (|3.37|) that 



(3.38) Ms)\\i,,^ds<Z + Ms)r,ds<Z. 

We now estimate ||9fp||o for £ = 0,1. By the equations (|3.35p - (|3.36p . the estimates (|3.33p . (|3.38p 
and (|3.5p . and the trace theorem, we find that for ^ = 0, 1, 



lia,V(.)|ii2(s^) + ||iafpi(.)|ii2(s^)ds 

(3-39) < 1^ \\didsus{s)\\l,^^^ + ||af5(s)||i2(s) + \\d'Ms)\\lds 

<Z+ f £{s)V{s)ds<Z. 
Jo 

By Poincare's inequality on il+ (Lemma IA.5p . ()3.33p and ()3.39p . we have 

r \\dfp+{s)\\l,^^^^ ds = f ||a,V(s)||i2(n^) + ||V9fp+(.)||i2(^^) ds 
(3.40) ''^ 



(3.41) 



< ||afp+(s)|ii2(s^) + ||v9,V(5)lli2(f,^)d.<z. 

On the other hand, by the trace theorem and (|3.39p ~ (j3.40p . we have 

ii5fp-(«)iii^(E_)^« < [ ii9fp+(«)iii2(s_) + mp{smi,^^_^ds 

ll^'p+(«)llii(c+) + III5'p(«)llli2(2_) ds < Z, 
and so again by Poincare's inequality on as well as (|3.33p and (|3.4ip . we have 

f \\dfp.{s)\\l,^^ ^ds= f \\dfp-{s)\\h^n ) + l|V5,V(s)||i2(f, )ds 

(3.42) 

< r ||a,V(^)lli2(s_) + I|v9fp_(s)|ii.(^_) ds < z. 



< 



In light of (|3.4Up and (|3.42p . we may improve the estimate (|3.33p to be 
(3.43) 



TE l|5^(^)lll-2m + E WdfpmUeds < Z. 



1=0 1=0 

Step 4 - Remaining V estimates for u and p 
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Next, we want to estimate \\pdfu{s)\\'^^j^-^^, and ||9|p||g. To do this, we shall use a weak formu- 
lation of the problem (jl.20p . It follows that 

(3.44) ^ 

- P+9idtV+, V3,+)+ + M aidh-^'fs)- for all ip G oH^in). 

Here (•,•)* denotes the dual pairing between (o-ff^(^))* ^.nd oH^i^), (■r)-i/2 denotes the dual 
pairing between //^/^(S) and H~^^'^{T,), {■,-)± is the inner product on I]±. By the estimates 
(133]) . (13311) and (|333]) . we see from the formulation dOD) that 



(3.45) f\\pdMs)\\lHiin)rds<Z. 

Since oH^{Q) C oH^{Q), the usual theory of Hilbert spaces provides a unique operator E : 
{oH^{n))* {oH^{n))* satisfying the property that Ef\^HUn) = f and that = 
ll/ll(o-H'^(r2))* fo'^ all / ^ (o-f^CT(^))*- Using this E, we regard the element of {oHl{Q))* as an element 
of {oH^{^))* in a natural way. Hence, we may also rewrite (j3.45p as 

(3.46) j\\pdMs)\\lHHn)rds<Z. 

With the estimate (|3.46|) in mind, since the pressure can be viewed as a Lagrange multiplier to 
(j3.44p , applying Lemma IA.91 we obtain 

(3.47) f\\dMs)\\Us<Z. 

Jo 

Now we derive the boundary estimates for dfu and dfp. First, by (j3.43p . we have 

(3.48) 1^ \\^*dMs)\\l^i/2^^) ds < 1^ \\dMs)\\lir2i^) ds < ||5^(s)||? ds < Z, 
and then by the incompressibility condition and (j3.48p . we get 

(3.49) 1^ 11535^3(^)11^-1/2(2) ds = 1^ \\d,d!ui{s) + 925^2(s)||^-i/2(s) ds < Z. 

For the term lud^d^Uil, i = 1,2, we use the first two identities of the dynamic boundary conditions 
on S_, along with the estimates (|3.5p and (|3.48p . to see that 



(3.50) 



\\ifid0tu,{s)}r^_,,,^^_^ds = yji5^i(s)ii^_v2(s„) + iiiM5N3(s)i(s)ii^-v2(s_)ds 

<Z+ f S{s)V{s)ds<Z, i = l,2. 
Jo 



Similarly, we have 

(3.51) 1^ l|9392n.,+(s)||^_,/,(^^) ds < Z. 
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Next, using the third identity of the dynamic boundary condition on E, along with the estimates 
(j3.5p and (j3.49p . we obtain the bound 

:i2„ /„M|2 , ||[ra2„/ Mi||2 



l^tP+(s)llH-i/2(s+) + III5iP(«)lllH-i/2(E_)^^ 

(3-52) < f ||a,V(s)ll-i/2 + I|93a,2n3(s)||^-i/2(s) + \\dh{s)\\l„ ds 



< z + 



f 8{s)V{s)ds<Z. 
Jo 



We now combine i^M), (l338|) . (lOSD . ^Ml^<^iMi) and (I3:i8]l - (l332]l . and then use the smallness 
of S to absorb the J term onto the left. This yields the inequality 

(3.53) fv{s)ds<Z, 

Jo 

which completes the dissipation part of the estimates. 
Step 5 - Conclusion 

Consequently, by ()3.28p . ()3.53p and the definition of Z, we have 

(3.54) £{t)+ I V{s)ds<Ce£{0) + Ce [ £{s)V{s)ds + eC I V{s)ds + Ce[ \\v{s)\\lds. 

Jo Jo Jo Jo 

By taking e and 6 sufficiently small, the previous inequality implies that 

(3.55) £{t)+ I V{s)ds <C£{0) + C [ \\r]{s)\\lds. 

Jo Jo 

This is (j3.23p and we conclude our proof. □ 

3.2. Energy estimates without surface tension. In this subsection, we will derive the nonlinear 
energy estimates for the system (|1.26p . For this, the definitions of the energy and dissipation rely 
on the linear energy identity of the homogeneous form of ()1.26p : 

According to this energy identity and the structure of the equations, as in [l^ 123] . we define the 
energies and dissipations as follows. For any integer > 3 we define the energy as 

2Af 2Af-l 

(3.57) £2N ■■= Yl (ll^'^llI^-2, + \\d{v\\lN^2j) + E \\dip\\lN-2,-i 

j=0 j=0 
and the corresponding dissipation as 

2Af 27V-1 
'D2N := Yl ll^'«ll47V-2,+l + E \\9tP\\lN-2j 

(3.58) 

+ ll'7ll4iV-l/2 + ll^t^ll4iV-l/2 + \\^tV\\lN-2j+5/2- 

We also define 

(3-59) T2N ■■= lhll4Ar+i/2- 

We will now derive our a priori estimates. We assume throughout this subsection that £2N{t) + 
J'2N{t) < 5^ for some sufficiently small 6 > and for all t in the interval in which the solution 
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is defined. We will implicitly allow 6 to be made smaller in each result, but we will reiterate the 
smallness of 6 in our main result. 

3.2.1. Energy evolution in geometric form. We can not use the equations of perturbed form (jl.26p to 
derive the energy evolution of the pure temporal derivatives because of lack of control of the highest 
temporal derivative of p. Instead, to control the temporal derivatives we employ the following 
equivalent geometric formulation. The geometric form is the equations directly transformed by 
composition from the original problem (jl.l2p : 



(3.60) 



(n,??) \t=o= (no,r?o)- 

Here we have written the differential operators Vyi,div_4 and with their actions given by 
(V_Af)i := Aijdjf, div^X := AijdjXi, and A^/ = div^ V^/ for appropriate / and X. We write 
Sj({p, u) := pi — fjMj^u for the stress tensor, where (JS)_Au)ij = AikdkU^ + AjkdkU^ is the symmetric 
^-gradient. Note that if we extend div_4 to act on symmetric tensors in the natural way, then 
div^ S_a(j), u) = — //A^n + Vj^p for those u satisfying div_4 u = 0. 

Applying the temporal differential operator df to (]3.60p . the resulting equations are the following 
system for {dfu,dfp,dfri), 

pdtidfu) - pWdsidfu) + pu ■ V^(afn) + dwASAidfp, dfu) = F^'^ in Q 



pdtu — pWdsu + pu 




- pAau + VaP = 


in Q 


divA u = 






in Q 


dtT] = u ■ Af 






on S 


Sa{p+,u+W+ = P+ 






on S+ 


M=0, lSAip,n) 


W- = 


Ipj gv-M^ 


on S_ 


u = 






on 



(3.61) 



where 



rn 



SA{d',p+,d',u+)M+ = p+gdir^+M+ + Ff 

Idlui = 0, ISAid'tP, dlu)iN. = M gdfv-^f- - f'J 

dfu- = 



in ft 
on S 
on 

on S_ 
on Sft, 



(3.62) 



f/'^= crpdrwdt"'d3Ui+ c^dt'^dtO^drKdaUi 

0<m<£ 0<m<e 

- Y cr{pdr{ujAjk)dt"'dku,+drAkdt"'dkP) 

o<m<e 

+ Y CrpdrAjkdt'^dkiAsdsUj+AjsdsUi) 
o<m<e 

+ Yl CrAjkdk{drAsdt"'dsUj + drAjsdt"'ds 



o<m<e 



+ dfdtO^KdsUi + Ajkdk{dtAisdsUj + dlAjsdsUi), 



(3.63) 



(3.64) 



F' 



- Y CrdrAijdt""djU^ - dfAijdj 

0<m<£ 

Y cr5rv7?+(af— r?+-af->+) 



o<m<e 



Y crp+ (dr{j\fj,+A,s,+)dt'^dsu,,^ + dr{j\fj,-,A,s,+)dt'^d, 



0<m<^ 
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o<m<e 

(3.65) - , 

0<m<l 



t 



U. 



(3.66) 

0<m<£ 

We present the estimates of these nonhnear perturbation terms p^'^ , F^'"^ , F^'^ and F^'^ in the 
following lemma. 

Lemma 3.5. Let the F^'*, i = 1, 2, 3, 4 he given as above, then for < £ < 2N , we have 
(3.67) WF'^X + mJF'^^)\\l + ||F^'3||2 + WF'^X < S.^V^n- 

Proof. Our perturbations F^'^, F^'"^, F^'^, and F^'^ have the same structure as those of [lOj. As 
such, the estimates in (|3.67p are recorded in Theorem 4.1 of |TOj. □ 

We now estimate the evolution of the pure temporal derivatives. 
Lemma 3.6. It holds that 

2N 



Y.\\^t<t)\\l+ f Y.\\^t<s)\\ids 

< ^2JV(0) + (£:2iv(t))'/' + f VS2N{s)V2N{s)ds+ f \\T]{s)\\lds. 

Jo Jo 



(3.68) ^=0 -"^ <?=o 



Proof. For £ = 0, 1, ... , 2N, we take the dot product of (|3.6ip ^ with Jdfu, integrate by parts over 
the domain and then integrate in time from to by using the other conditions in (j3.6ip and 
some easy geometric identities, we obtain the following energy identity: 

i / pj\dfu{t)\' + ^ f j /.J|D^9^|2 

(3.69) =\j pAcfM^t + f I Jidfu . F' + dfpF') - f I dlu ■ F' 

^ Jn Jo Jn Jo Jt. 



First, we estimate the left hand side of (j3.69|) . For this we write 

(3.70) J|D^afnp = \Bdfu\^ + (J - l)\Bdfu\^ + JiB^dfu + Odfu) : (B^^fn - B^fu). 
For the last term in the above, since 

(3.71) B_Ad^^u ± D9fn = {A^k ± 5^k)dkdiu^ + {Ajk ± 5jk)dkdiu\ 
we have that (using Lemma lA.ll to estimate ^ it /) 

(3.72) / J(B^afn + mlu) : (B^afn - B^fu) < v^^27V- 
Jq. 

On the other hand, we similarly have the estimate 

(3.73) f |J-l||Bafn|2<7^P2Jvand / p\J - l\\dluWt) < {82Nf''' ■ 
Jq Jn 
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This allows us to estimate the left hand side as 



(3.74) 



> \ I P\dtu{t)\^ + \ f I ^ipdlu\''-C{S2N{t)f'^-C f ^/£M^V2N{s)ds. 
^ Jn ^ JO JQ Jo 

We now estimate the right hand side of (|3.69p . For the term, by (j3.67p . we have 

(3.75) / / Jdfu-F^< [ \\dfu\\o\\J\\L^\\F^\\o < [ ^/&^V2n. 

Jo Jn Jo Jo 

For the F^ term, by again (j3.67p and the trace theorem, we have 

(3-76) - /* / dfu • F3 < f \\dfu\\ms)\\F% < f \\dlu\\^^/£2NV2N < f V£^V2N. 

Jo Jy. Jo Jo Jo 

Next, we estimate the F^ term. Notice that we can not control d^^p by I?2Af; hence we need to 
integrate by parts in time to find 

(3.77) f [ JdfpF' = - f [ dt'pdtiJF^)+ [ (dt'pJF'm- [ {dfpJF'm. 
Jo Jn Jo Jn Jn Jn 

Then by (|3.67p . we may estimate 

(3.78) - r / dt'pdt{JF^)< ['\\dt'p\\o\\dt{JF^)\\o< f V^V2 

Jo Jn Jo Jo 

Also, it is easy to deduce that 

(3.79) / {dt'pJF'm - [ {dt'pJF^m < £2NiO) + {£2N{t)f''. 
Jn Jn 

Hence, we have 
(3.80) 

Jo Jn 

For the last line of (j3.69p . by the trace theorem and Cauchy's inequality, we have 



JdlpF^ < £2n{0) + {£2N{t) f/^ + /* V£^'D2N. 

Jo 



(3.81) 



Jo J S-f Jo 



<C I \\dtri{s)Mdtu{s)\\ods<C, 1^ \\dfr,{s)\\lds + el^\\dfu{s)\\lds. 



Using Korn's inequahty together with the estimates (|3.74p . (|3.75p . (|3.76p . (|3.80p . and (|3.8ip . we 
find that we may take e sufficiently small to deduce from ()3.69p that 



WMtm+ / WMs)\\ids 

(3.82) ^ ^ 

<£2N{0) + {£2N{t)f^+[ ^/£Ms)V2N{s)ds+ f ?7(s) Hg ds. 

Jo Jo 



Now taking ^ = in (j3.82p . we have 
\\nml+ f\Hs)\\\ds 



(3.83) 



£2N{^) + {£2N{t)f'^ + f ^£2N{s)V2N{s)ds+ \\ri{s)\\l ds. 

Jo Jo 
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For a. = 1, . . . , 2N the kinematic boundary condition (jS-GOP g and the estimate (j3.67p imply that 

(3.84) liafr/llo < ||af-in3||i2(s) + \\dl-\u ■ < Wdf^Wi + V^V^- 

Plugging ([3:81) into iK82\i . we obtain 



(3.85) 

<£2N{0) + {£2N{t) f^+l V£M^V2N{s)ds+ I \\dt' u{s)\\l ds . 



Consequently, chaining ()3.85p and ()3.83p . we get ()3.68p . □ 

3.2.2. Energy evolution in perturbed form. We shall now use the perturbed equation (jl.26p to derive 
the energy evolution of the mixed horizontal space-time derivatives. To proceed with, we present 
the estimates of defined by (|1.27p - (jl.3ip in the following lemma. 

Lemma 3.7. There exists a 9 > so that 

(3.86) WV^'^-'GX + IIV^^-^G^II^ + ||Vr-'G3||?/2 + IN'f-'CX/, < S^^', 

||f^4iV-2/-Yl||2 I \\fj4N-2r<2\\2 , ||f^4Af-2^3 1|2 , ii f^4Ar-2^4 1|2 
II V ^ llo + ll^ ^ llo + ll^* '-^ lll/2 + ll^* lll/2 

(3.87) + WV^^'-^dtGX + WV^'^^'dtGYo + l|Vr~'5tG3||2/2 + WVt'^'^GYy, 

< £InV2Ni 



(3.88) II V G llo + IIV G ||o + ||V^ G II1/2 + 11^* G II1/2 ^ '^2Ar'^2Af + i2Af>2Ar- 

Proof. Since our perturbations G^,G^,G'^, G^ have the same structure as those of [10], the estimates 
dMSIl-dSSH]) follow from Theorem 3.2 of [lOj. In the statement of Theorem 3.2 of [10], the left 
hand side of ()3.88p is bounded by £2n'^2N + ^^2N with 

2 

(3.89) /C := ||Vn||ioo + ||V\||ioo + ^ \Duif^2(j^y 

1=1 

However, in the present case we may use the Sobolev embeddings to estimate fC < £2N, and so we 
have dSSH]). □ 

Now we estimate the energy evolution of the mixed horizontal space-time derivatives. 

Lemma 3.8. There exists a 9 > such that for any k > 0, there exists a constant > so that 



(3.90) 1"!-''^ I"!-"^^ 



< 



£2n{0) + ^ {£In'^2n{s) + ^/£^:F2n{s) + kF2n{s) ds^ +C^J^ 



MsWods 



Proof. Since the boundaries of are flat we are free to take the time derivatives and horizontal 
derivatives in the equations ()1.26p . Let a G N^"*"^ so that oq < 2N — 1 and |a| < AN. We apply 
d" to ()1.26p ^. multiply the resulting equations by d'^u, and then integrate over using the other 
conditions in ()1.26p : using trace estimates on the resulting equality, we find that 



2dtJ^^^ ' 2 J^^^ ' 

3.91 



a„,i2_ / • a"G^ + a>5"G2 - / d'^u-d'^G^ 
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We will estimate the right hand side of (j3.9ip . First, we estimate the G^,G^, terms. Assmne 
that \a\ < AN — 1, then by the estimates (|3.87p - (|3.88p in Lemma 13.71 we have 



(3.92) 



/ • d'^G^ + ayd'^G^ < ||9"n||o||a°G^||o + ||5>||o||a"G^||o 
Jn 

Again by (j3.87p - (j3.88p . together with the trace theorem, we have 

(3.93) - j^'^u-d'^G^ < ||9%||i2(2)||9°G3||o < V^^J£'^J^V2N + £2nJ'2N. 

Now we assume that |a| = 4N. Since oq < 2N — 1, we have ai + Q2 > 2, so we can integrate 
by parts in the horizontal directions. We write 9" = d^d'^ so that I7I = 4N — 1 and |/3| = 1. 
Integrating by parts and using the estimates (j3.87p - (j3.88p . we find that 

[ a"n-a°G^ = - I d'^+^u-d'^G^ <\\d'^+^u\\o\\d'^G^\\o 
Jn Jn 

^ \/T^2N\J £2n'^2N + £2N^2N- 

For the G^ term, we do not need to integrate by parts: 

(3.95) < WdyWoWG'^Wi < V^\/£2N^2N + £2nJ'2 

For the G^ term, we use the trace theorem to estimate 



(3.94) 



N- 



S^+^-u • d^G 



<W^^u\\h-V.^^)WGX/2 



(^•96) < Wu\\h.,2^^)WGX/2 < \\d"uUd''G%/2 



^ \/T^2N\J £2N^2N + £2N^2N- 

Now we turn to estimating the last two terms in ()3.9ip . By the trace theorem and Cauchy's 
inequality, since uq < 2N — 1 and |q;| < AN, we have 



(3.97) 



p+gdy+d^us,+ + ll^ Mgd''V-d''us<G\\dy\\_y2\\d"u3\\Hir. 



<c,iiarMiL-2.o-i/2+^ii5"^ii?- 



Consequently, in light of (j3.92p - ()3.97p . we may integrate (|3.9ip from to t and sum over such 
a. Using Korn's inequality and taking e sufficiently small, we find that 



(3.98) 



\a\<4N \a\<4N 
ao<2N ao<2N 



+ / {£2N^2n{s) + ^/V2n£2nJ'2n{s)) ds + E W^?' r]{s)\\ 

If ao = 0, we use Sobolev interpolation to bound 

(3.99) Pr^ll4Af-2ao-l/2 = ll^ll4Af-l/2 < C'^^lkllo + '^11^11^+1/2 = C'^ll^llo + l^^2N ■ 
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If 1 < ao ^ 2A^ — 1, we use the kinematic boundary condition, the trace theorem and ()3.87p . to 
obtain 

ll^r'''?ll4Af-2ao-l/2 — "^^^3||^4JV-2qo-1/2 + \\di° ^G^||^4]V-2ao-l/2 

(3-100) <\\dT ^^ll4iV-2(ao-l)-2 + ll^r" ^G'^ll^4]V-2(c«o-l)-5/2 

<l|9r-^n||2^_2(,„_i)_2+^:2V2iV. 

Hence, by (|3:99]) ~ (l3lTO . we deduce from (l3:98D that 

|a|<4iV *^ |a|<4Af 

ao<2A' ao<2A' 

<£2n{0)+ / £2i^V2N + V£^J'2N + ^iT2N ds 

Jo 

+ Cj' Ms)\\lds+ I Il^'^(^) 114^-2^-2^^5- 

Jo Jo o<£<2Af-2 

To conclude, we combine (jS.lOip and (j3.68|) of Lemma 13.61 renaming 9; This yields (j3.90p . □ 

3.2.3. Transport estimates. Now we derive estimates for rj by using the kinematic transport equa- 
tion 

(3.102) dtr] + u-V^r] = U3 in S, 

where u • V*r/ = uidirj + U2d2r]. 

Lemma 3.9. For any e > 0, there exists a constant > such that 

2N 



(3.103) 



Y \\9fvit)\\l_2e + ^2N{t) < f27v(0) + J-2iv(0) + / V£^iV2N + ^2Jv) ds 

Jo 

+ e f{£2N + J'2N)ds + Ce f V Wu{s)\\lds. 
Jo Jo I 



|a|<4Ar 



Proof. Let a G Applying 9° with |a| < AN to (|3.102p . we obtain 

(3.104) dtid'^rj) + u • V*(9"r7) = a^ug + R"- 

Here the remainder term i?" is a sum of terms d^u ■ V*(3"~^r/) with 7^ /3 < a. As in Lemma [3. 7^ 
we may bound this term as 

(3.105) ll^"lli2(S) ^ ^2W(^27V +-^27V). 

We multiply (j3.104p by d°'r} and then integrate over S; by the Sobolev embedding theory on S, the 
trace theorem and (j3.105p . we obtain the bound 

1 d 



2 dt 



(3.106) 

< ||V,n||ioo(2)||9"ry||2 + (||5°n3||i2(2) + ||i?"||L2(s))||5-77llo 

< hMd'^vWl + (P^^slli + ||i?°||L2(E))l|5"r/||o 

< V£^{'^2N+J^2n) + Vs^Wd'^uWi. 
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(3.108) 



Integrating the previous inequality in time from to t, summing over a, we find that 

u{s)\\i ds. 

e=0 •'^ \ \a\<4:N ) 

Now we define the operator J = ^/T^^'Kl. We apply J^^+^I"^ to ()3.102p . multiply the resulting 
equation by i7^^'^^/^r/, and then integrate over (0, t) x S; using the commutator estimate of Lemma 
[OOl with s = 4iV + 1/2, we find that 

= \ |^(5ini + 92^2)1:^^^+^/%!' + £ {j'^'^^'^u, - [j^'^^'l^u] . V.r?) J^^^^l\ 

< \\V.u\\L^ij:)\\j'''+'/'v\\l + (|I^'^+'/'^3||l^(e) + l|V.n||^o.(s)||J^^-i/2V.r?||o 

^ lkl|4|hll4Ar+l/2 + (ll^^^^^3||l + II^^IU ||??||4Af+l/2 + ll^^^^^3 1| 1 ll^ylls) ||??||4Af+l/2 

< V£^{V2N + J^2n) + ^/^\\^''u\\l. 

Integrating the above in time from to t, we find that 

(3.109) T2N{t)<T2N{0)+ f ^/£^{V2N+J'2N) + ^^\\d''u{s)\\ids. 

Jo 

Summing (j3.107p and (j3.109p and applying Cauchy's inequality then yields ()3.103p . □ 

3.2.4. Full energy estimates. We will improve the estimates derived previously. To simplify nota- 
tion, we define the "horizontal" energy as 



2N 

2 

21 



(3.110) £2n:= E l|5"^llo + Ell^Mll 

|a|<4Ar £=0 

and the horizontal dissipation as 



(3.111) V2N-= 11^" 

|a|<4Af 



nil?. 



First, chaining ()3.90p and ()3.103p . we obtain that for any e > 0, there exist a constant > 
and a universal Cq > such that 

£2N{t)+J^2N{t)+ [ V2Nds 

Jo 

(3.112) < Cs{£2n{0) + J'2n{0)) + Ce [ £^n1^2N ds + Ce [ ^/£^T2N ds 

Jo Jo 

+ eCo / {£2N + J^2n) + CeKT2N ds + C^Ce / \\r]{s)fods. 
Jo Jo 
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If we set K = (eCo)/Ce, then we may rewrite the inequahty above as 

£2N{t)+J^2N{t)+ [ V2Nds 

Jo 

(3.113) < Cs{£2n{0) + -F2iv(0)) + Cs f £^nV2n ds + Ce f y/£^F2N ds 

Jo Jo 

+ 2eCo [ {£2N + J'2N)ds + Ce [ Ms)\\lds. 
Jo Jo 

To conclude our estimates, we shall replace the left hand side of (j3.113p with the full energy and 

dissipation. Hence, in what follows we want to show that £2n is comparable to £2n and that 'D2N 

is comparable to 'D2N, provided that 5 is sufficiently small. We begin with the energy estimate. 

Lemma 3.10. It holds that 

(3.114) £2N < £2N. 

Proof. We compactly write 

n-l 

(3.115) Z2N = ^ \\diG^\\lN-2j-2 + ll'^G'^ll4Ar-2i-l + ll^<^^ll4Af-2i-3/2- 

j=0 

Note that by the definitions of £2N, we have 

2N 

(3.116) Wd^uWl + \\dN\lj,.2, < £2N. 

3=0 

Now we let j = 0, . . . , 2N — 1 and then apply to the equations in (|1.26p to find 
-I^Adiu + Vd{p = -p4^^u + dlG^ in 



(3.117) 



dWdlu = diG^ inn 



ldi,p+I - //+D(a^u+)]]e3 = p+gdlv+es + d^Gl on S+ 

ldiu\ = 0, IdlpI - pB{diu)\e3 = M gdlrj^e^ - djCl on S_ 
dfU- = on Tif). 

Applying the elliptic estimates in Theorem IA.7I with r = iN — 2j > 2 to the problem ()3.117p and 
using ()3.115p - (|3.116p . we obtain 

\\dlu\\l^-2j + \\dlp\\lN-23-l < ll^'+'^llL-2,-2 + \\01G'\\In^2j-2 
(3-118) + ||c^G^|l4Ar-2j-l + ll'5^^ll4Af-2j-3/2 + ll^t ll4Ar-2j-3/2 

^ ll^^^^ll4Af-2{i+l) +£2N+2:2N- 

We claim that 

(3.119) £2N<£2N+Z2N- 

We prove the claim (|3.119p by a finite induction based on the estimate (|3.118p . For j = 2N — 1, 
we obtain from (|3.118p and (|3.116p that 

(3.120) Wd^'Ml + Wd^-'pfi < Wdr^fo + £2N + Z2N < B2N + Z2N. 
Now suppose that the following holds for 1 < i < 2N 

(3.121) \\dr-'u\\l, + \\dr-'p\\l,., < £2N + Z2N. 

We apply ()3.118p with j = 2N — (£ + 1) and use the induction hypothesis ()3.12ip to find 



(3.122) \\dr'^'^'\\\\,^,^ + < \\dr-'uf2, + £2N + Z2N < £2N + 



2N- 
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Hence by finite induction, the bound (j3.12ip holds for all I = l,...,n. Summing (j3.12ip over 
i = 1 , . . . , n and changing the index, we then have 

2N-1 

(3.123) Yl 11^' 

j=0 

We then conclude the claim (j3.119p by summing (j3.116p and (j3.123p . 

Finally, using ^M\\ in Lemma EZl to bound Z2N < ('?2Af)^+^ we obtain (|3.114p if we take 
£2N < with 6 small enough to absorb onto the left. □ 

Now we consider a similar estimate for the dissipation. 
Lemma 3.11. It holds that 

(3.124) V2N <'D2n + £2nJ^2N. 
Proof. We compactly write 

(3.125) Z2N = ||V^^-^Gi||g + ||V^^-^G2||2 + ||vr"'G3||2/2 + ||Vr-'G^II?/2 + l|V^^^ 
First, by the definition of 'D2N and Korn's inequality, we obtain 

(3.126) llVr^ll? <^2iV. 

Notice that we have not yet derived at estimate of rj in terms of the dissipation, so we can not 
apply the two-phase elliptic estimates of Lemma [A. 71 as in Lemma [XTUJ It is crucial to observe that 
from (|3.126p we can get higher regularity estimates of u on the boundaries S = S+ U S_. Indeed, 
since Ti± are flat, we may use the definition of the Sobolev norm on and the trace theorem to 
see from (I3.126P that 

||^«||^4iV-2j + l/2(s) ^ l|t^^llL2(2) + II ^■'9^n||^i/2(j^) 

<mu\\l + \\Vt''-''d{u\\i<'D2N. 

This motivates us to use the one-phase elliptic estimates of Lemma IA.81 
Let j = 0, . . . , 2N — 1, and observe that {dl.u+,dlp+) solve the problem 

-fi+Adiu+ + Vdip+ = -p+di^\+ + d{Gl in n+ 
(3.128) { divd{u+ = diGl in n+ 

= dju^ on S, 



and {dlu-,d^p-) solve the problem 
(3.129) 



-fi-Ad{u- + Vd{p^ = -p-d{^\- + did in 
div dfu^ = dfG"^ in 



c^u_ = on S_ 

= on Tib. 

We apply Lemma IA.8I with r = iN — 2j-\-l to the problem ()3.128p for , p+ and to the problem 
(13.1291) for n_, p_, respectively; using (I3.125p . (I3.126p . (I3.127P and summing up, we find 

\\diu\\lr,.2j+i + l|Vc^p||L_2j-i 
(3.130) < ||a^^n||l^_2,_i + \\diG'\\l^_2,_, + \\diGYm-2j + l|5^'«ll^4iv-2,+i/2(s) 

^ \\di^^u\\lj^^2j-l + Z2N + V2N- 
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We now claim that 

2N 2N-1 

Z2N + T^2N- 

3=0 i=o 

We prove the claim (j3.131|) by a finite induction as in Lemma I3.1UI For j = 2N — 1 , by (|3.125p and 
(|3.13U|) . we obtain 

(3.132) Wd^'Ml + II Va^^- VII? < ll^r^^ll? + + f)2N < Z2N + f)2N- 

Now suppose that the following holds for 1 < i < 2N: 

(3.133) \\0r~V2i+l + l|Va2^-Vll2^-l < Z2N + t>2N. 

We apply (j3.130p with j = 2N — (^ + 1) and use the induction hypothesis (j3.133p to find 

^3^3^^ ll^t ^Il2(m)+1 + ll^^t Pll2(m)-1 ^ ll^t U\\2i+1 + 22N+T)2N 

< Z2N + T^2N- 

Hence the bound (j3.133p holds for all / = 1, . . . , 2N. We then conclude the claim (|3.13ip by summing 
this over £ = 1, . . . , n, adding (|3.126p and changing the index. 

Now that we have obtained (|3.13ip . we estimate the remaining parts in T)2n- We will turn to 
the boundary conditions in (jl.26p . First we derive estimates for 77. For the term c^?] for j > 2 we 
use the boundary condition 

(3.135) atr/ = na + on S. 

Indeed, for j = 2, . . . , 2iV + 1 we apply to (I3J35D to see, by ([TOT]) and ([3l^ . that 

ll'^'?ll47V-2i+5/2 ^ ll^t ""311^^-2^+5/2(2) + ll^t G'^ll4Af-2j+5/2 

^ ll^r «3|l2n^2(i-l)+l + 11^" G^ll4iV-2(j-l)+l/2 ^ Z2N+V2N- 

For the term dtrj, we again use p.l35p . (|3.13ip and (|3.125p to find 

^^^^ ll^t^ll4Ar-l/2 ~ ll'"3||^4JV-i/2(s) + 11^ ^^"^ ll4Ar-l/2 

^ ll'"3|l4Af + l|G'^ll4Ar-l/2 ^ ^2Af + ^2Ar- 

For the remaining ry term, i.e. those without temporal derivatives, we use the boundary conditions 

(3.138) pgri+ = p+ - n+d3U3,+ - on S+ 
and 

(3.139) Ipj gri^ = M " Im^s^sI + G|_ on S_. 

Notice that at this point we do not have any bound on p on the boundary E, but we have bounded 
Vp in 0. Applying V^, to p.l38p and (|3.139p . respectively, by (|3.13ip and (|3.125p . we obtain 

^^^^ l|V*^ll4Af-3/2 ^ l|V*Pll^4iV-3/2(s) + l|V*t^3'U3||^2n-3/2(s) + ll^*G'^ll2n-3/2 

< l|Vp||L-l + ll^^3|ll7V+l + l|G^ll4JV-l/2 ^ ^2N + V2N ■ 

Since /rp2 = 0, we may then use Poincare's inequality on Ti± to obtain from (j3.140p that 

(3-141) ll^ll4Af-l/2 ^ ll^llo + l|V*?/||4iv_3/2 < II V*?7||4jv_3/2 ^ ^2Af + ^2Af- 

Summing ()3.136p . ()3.137p and ()3.14ip . we complete the estimates for r/: 

n+l 

(3-142) ll^ll4Af-l/2 + ll^t^ll4Af-l/2 + ^ ll^^ll4Ar-2j+5/2 ^ Z.2N + 'D2N ■ 

j=2 
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It remains to bound \\d{p\\o. Applying j = 0, . . . , 2iV - 1 to (IXT38]l - (l3J39D . by (imi]) . 
(I3T12D and (f3l^ . we find 



(3.143) ... 

< \\&lrj\\l + WdlusWl + WdiGYo < Z2N + V2N. 

By Poincare's inequality on 0+ (Lemma IA.5P and (I3.13ip and (]3.143p . we have 

(3 144) W^tPWW,) = \\dip\\l.^n,) + W^^tPWhin,) < \\&ip\\W,) + l|V5^p||i.(o,) 

< Z2N + '^2N- 

On the other hand, by the trace theorem and (|3.143|) - (|3.144p . we have 

< Z2N + T^n, 

SO again by Poincare's inequality on as well as (|3.13ip and (|3.143p . we have 

(3 146) \\^tP\\mi^_) = WdipWl^n^) + W^mUn^) < + l|Va^p||i.(f,_) 

< Z2N + T^2N- 

In light of (|3.144p and (|3.146p . we may improve the estimate (|3.13ip to be 

2N 2N-1 

(3.147) ^\\di 

\\9{p\\lN-2j<22N+V2N- 

j=0 j=0 

Consequently, summing p.l42p and (j3.147p . we obtain 

(3.148) Vn<Z2N+'D2N. 

Thus, using dSSZD-dMHl) in Theorem EZl to bound Z2N < {£2nYT>2n +£2nJ'2N , we obtain (i3T2i]l 
provided that £2n < with 6 small enough to absorb onto the left. □ 

Now we are ready to state the full energy estimates. 

Theorem 3.12. Let A* be defined by dOH]) . and let (n,p,7/) solve (fLMIl . If £2N {t) + T2N (t) < 6^ 

for sufficiently small 6 for all t € [0, T], then we have that 

(3.149) 

£2N{t)+J^2N{t)+ [ V2N{s)ds 

Jo 

<C{£2N{0)+J'2Nm + ^ [ {£2N+J'2N)ds + C [ y^T2N ds + C [ Ms)\\l ds 

^ Jo Jo Jo 

for all t G [0,r]. 

Proof. By (I3.114P and (j3.124p . we deduce from (13.113P that for any e > 0, there exist > and 
Ci > such that 

£2N{t) + J'2N{t) + f V2N{s)ds 

Jo 

(3.150) < Ce{£2N{^) + -F2JV(0)) + Ce f £^n'^2N ds + Ce [ ^/£^T2N ds 

Jo Jo 

+ eCi [ {£2N + J'2N)ds + Ce [ Hs)\\Us. 
Jo Jo 



38 



YANJIN WANG AND IAN TICE 



We fix the value of e by e = A*/(2Ci) with A* defined by (j2.33p . Taking 5 sufficiently small, we 
deduce from (j3.150p that there exists a constant C = Ca, so that 

£2N{t)+J^2N{t)+ [ V2N{s)ds 

(3.151) t t t 

<C{£2N{0)+T2Nm + ^ [ i£2N + J'2N)ds + C [ y/£^F2N + C f Ms)\\lds. 

^ Jo Jo Jo 

This is exactly ()3.149p and the proof is completed. □ 

4. Nonlinear instability 

4.1. Unified operator form. Note that the nonlinear problems (jl.20p and (|1.26p can be written 
uniformly as 



(4.1) 



pdtu — /uAti + Vp = in f2 

div u = div F'^ in Q 

dtf] — Us = — -63 on S 

{P+I - fJ'+'^{u+))e3 - {p+gri+ - cr+A*7?+)e3 = F| on 
M = 0, b^-//B(7x)le3-(M5??-+<^-A*r?-)e3 = -F3 on S_ 

n_ = on Sfe 
(u,r?) |t=o= {uQ,r]Q). 



Here, the nonlinear terms in the right hand sides of (j4.ip are given as follows: for (t± > 0, ()1.20p 
corresponds to ()4.ip with F^ = f, F"^ = and F^ = g, where / and g are defined by ()1.2ip - (|1.25p : 
for CT± = 0, (fT:26]l corresponds to gl|) with F^ = G^,F^ = {5ji - JAji)uj, i = 1, 2, 3 and = G^, 
where and are defined by (fL27l> . (fOOl) and (fOOl) . Note that for a± = 0, (jEJg is derived 
as follows: 

(4.2) div^n = <;=^ JAijdjUi = <;=^ dj{JAijUi) = <;=^ divu = divF^. 

This implies that divF^ = G^, with G^ defined by ([L28]). We can also check, using the definitions 
in (jl.lSp . that F^ ■ = —G'^ with defined by (|1.3ip as follows: on the boundary S, 

(4.3) f2 . 63 = {5ji - JAji)uj5i3 = u3-u-M = -G\ 

In order to use the linear theory we have developed in Section [2l we shall write the system 
()4.ip in the operator form as in [H El [3]. First, we want to reduce the divergence of velocity to 
be zero. This can be easily done by setting u = u — F"^. Then {u,r],p) satisfy (14. ip . replacing F^ 
by = - pdtF'^ + )uAF2, F2 by 0, and F^ by F^ = F^ + ijM{F^). Second, we reduce the 
tangential components of the stress tensor F?, i = 1, 2 to be zero. We choose the vector z satisfying 
the conditions (see Page 6 of [3], for instance, for details of how to construct z) 



(4.4) 



2+ = 0, 93Z+ = 0, — = (Ff +' ^+ vanishes near S_ 

z_ = 0, d^Z- = 0, —fi-d'^Z- = (— FI -^f 0) on z_ vanishes near E;, 

Then w = V x z satisfies 



w = 0, -p+{d3Wi + diWs) = F?+, dsws = on S+, 

(4.5) { w = 0, -Md3w' + diw^)j = -Ff_, dsws = on S_ 

divw = in Q, and w = on Sf,, 

and we have the estimate 

(4.6) lhll.< l|i^'llr-3/2, ''>2. 



Therefore, let u 
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pdfU — fj.Au + Vp = F^ := F^ — pdtw + pAw in Q 
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div u = 






in 


dtV -U3 = 






on S 


Stan{u+) = 






on S+ 


P+ - "^l^+d^u^,- 


h - P+9^+ + f^+A*r/+ 


= PL 


on S+ 


M = IStanm 


= 


on S_ 


\p - 2^53n3]] - 






on S_ 


U- = Q 






on Eft 



(4.7) 



(n,??) |t=o= (uo,??o), 

where we have written Stan{u) for the tangential part of S{p,u). 

We want to introduce a projection on divergence- free vectors to remove the pressure as an 
unknown. With fixed boundaries the pressure term has no influence on the projected equation; 
in the present case, we shall see that it is replaced by a gradient term which is determined by 
the other unknowns. Let V be the orthogonal projection in L^(r2) onto the space orthogonal to 
{V(j) I G °H^{n)}, where ^H^{n) := {u G H^{n) \ u\t:+ = 0}. Following Lemma 3.1 of [IJ, 
we can prove that V is a bounded operator on H^{Q.) for s > and that if iji) E H^^^{Q.), then 
V{'V(j)) = Vvr, where vr solves the problem 

(4.8) Avr = in vr = on S+, {irj = {(pj and {dsTrj = on d^n = on S^. 

We may regard vr as a harmonic representation of cp and we denote it as tt = ?^(0|s+, Il0l|s_) to 
illustrate its dependence on the boundary values of (p. 

Now we apply V to the momentum equation (j4.7P ;^ to obtain 



(4.9) pdtu + Au + Vnip+g7]+ - a+A,r]+, {pj c/r?_ + A,r/_) = VF^ - VUiF 
where A is defined by 

(4.10) kv = -pVAv + Vn{2p+d3V:i\j^^,2lpd3V3i\^ 
The domain condition of A is given by 

(4.11) div?; = in Vt, Stan{v) = on S_ 
Now we define the linear operator C by 

''u\ ( p~^ (An + VHip+gr]- 



Ivj = lStan{v)j = on S_, u = on Sfe 



-U3 



+ c 



(4.12) C 

which then allows us to rewrite (j4.7p as 

(4.13) dt 

where the nonlinear term ^ is given by 

(4.14) m := p-^ [VF' - Vn{F. 
Employing Duhamel's principle, we can then solve (j4.13p via 

(4.15) 



(T_A^,7?_)) 



u{t) 

We can then go back to u = u + w + F"^ to see that 







ds. 



(4.16) 



+ 



ai(s) 




ds. 
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4.2. Unified estimates. We now define the norm |||-|||q appearing in the Theorems 1 1 . 1 1 and 1 1 . 2 1 by 

I I ii\ I I \ 

(4.17) 







fIIo + ll^l 



0- 



We define the term 
(4.18) 



loo appearing in the Theorem 11.11 by 

:= VS, 



00 



where £ is defined by (j3.2p . We define the term 



II 00 



appearing in the Theorem 11.21 by 



(4.19) 



:= \/ £2N + J^2N 



00 



for an integer > 3, where £2N and T2N are given by (|3.57p and (|3.59p . respectively. We also 
define another norm ||-|| by 



(4.20) 



Ml + Mil 



Remark 4.1. Note that the expression for 



also involves sums of norms of p and its 



00 



II 00 



as 



temporal derivatives, and as such is not actually a norm for u and rj. If we were to view 
acting on triples {u,p,r]) in the obvious way, then this quantity would function as an actual norm. 
Here we abuse notation in this manner in order to suppress the appearance of p, which we want to 
view as being determined by u and rj in a nonlinear fashion. The reader troubled by this abuse of 



notation could simply replace the notation 
the same manner. 



by a nonlinear functional <B{u, rf) defined in 



00 



We now restate the main results of Sections [2] and [3] in our new notation. Note that our notation 
has allowed us to unify the cases (T± = and a-i- > 0. 

Proposition 4.2. Suppose that {u,p,r]) is the solution to (ll.20p with surface tension or is the 



solution to (|1.26p without surface tension. Let the terms \1-\Iq, 
()4.20p . Then we have the following. 

(1) Let A be defined by ()2.32p . Then we have 

(n{0)\ ' 



II 00 



and 



be defined by (gH])- 



(4.21) 



< Ce 



(2) There is a growing mode 



(4.22) 



satisfying 



2At 

V* 

.At / u- 



(u{0)\ 



< 00, and 



00 



(3) There exists a small constant 6 such that if 



uit) 
r]{t) 



< 6 for all t € [0,T], then there 



00 



exists Cs > so that the following energy estimate holds for t G [0, T] : 



(4.23) 



u{t)\ 

v{t)J 



<Cs 



00 



(u{0) 



+ Cs 



2 

+ 

00 
u{t)\ 

v{t)J 







uit) 

7]{t) 



2 

3 

ds + Cs 

00 



00 



ds 

uit) 
r]{t) 



ds. 



Here, in the statements of 2 and ?>, \ > Q is a number satisfying ^ < A < A. 
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Proof. Statement 1 follows from Theorem 12. 7[ Statement 2 follows from Theorem 12. 4i Statement 
3 follows from Theorem 13.41 by taking A = A for the case with surface tension and from Theorem 
13. 121 by taking A = A* for the case without surface tension. □ 

Now we will state an estimate for the solution to (j4.ip . 

Proposition 4.3. Suppose that (u,p,rj) solve (14. ip and that u = u — w — F"^ with w constructed 
in Section \4.1\ Then we have the following estimates: for both a± = and a± > 0, there exists a 

'u{t)\ 



small constant 6 such that if 



(4.24) 



00 

< C 



< 6 for all t £ [0,T], then 

2 



u{t) 

vit) 



+ C 



00 



.A{t~s) 



u[s^ 
vis] 



ds 



00 



for all t G [0,T]. 

Proof. From the formula (I4.16P and (14.2ip . we haye 

jc 



(4.25) 



uit) 
r]{t) 



(m 

U(0) 



< 



< 



w{t) + F^{t) 




w{t)\\o + \\F' 



+ 







ds 



+ 



=A(i-s) 



\ms)hd.s. 



On the other hand, by the construction of ^ and w in Section [4.H we can estimate 



(4.26) 






\\w\\ 


0<lkl|2<||i^'||l/2<l|i^'lll/2 + l|i^'l|2. 


and 












ll^lb < 


|F 


I2 + 


\\vn{Fl+,-Fl_)\\2 




< 


|F 


I2 + 


||5tU;||2 + ||'U;||4 + ||i^^||5/2 


(4.27) 


< 


\F' 


I2 + 


\\dtF% + \\F% + \\dtw\\2 + \\W\U + ||F'||5/2 + ll^(i^')ll5/2 




< 


\F' 


I2 + 


\\dtF% + Ili^'lU + ||i^'||5/2 + \\dtF%/2 + ||F'||5/2 




< 




I2 + 


\\dtF 2 + \\F 4 + \\F 5/2- 



Plugging g26])-(I12Zl) into K25\i . we obtain 
(4.28) 



v{t)J 



< \\F\t)h + 11^=^(^)111/2 + / e^(*-)(||Fi(s)||2 + \\dtF\s)h + \\F\s)\U + \\F\s)h/2) ds. 

Jo 

In both cases a± = and a± > 0, we may estimate the nonlinear terms as we did in Section [3l 
Indeed, we haye 

2 



(4.29) 

and 

(4.30) 



\F' 



1/2 ^ 



u{t) 
r]{t) 



00 



\F\s)\\2 + \\dtF\s)\\2 + \\F\s)\U + \\F'{s)h/2 < 



uis) 
r]{s) 



2 

00 



where we recall that |||-|||oo is given by (I4.18P when a± > and (14.19P when a± > 0. 
Substituting (lilMjl - flCTD into K2E\i . we get K24\i . 



□ 
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Thus far we have not elaborated on the local well-posedness theory for our problem that we 
developed in [23]. We now want to state a version of it specially tailored to our uniform operator 
setting. In the result we will refer to the "necessary compatibility conditions" required for the local 
well-posedness in our energy spaces. These are cumbersome to write out explicitly. In the case 
without surface tension we refer to Theorem 2.1 of [23] for their explicit statement. In the case 
with surface tension we refer to Theorem 2.5 of |23] for a lower-regularity version of the conditions; 
the conditions in the present case can be deduced similarly as those stated in [23j. 

Theorem 4.4. Suppose that the initial data {uo^rjo) satisfying the necessary compatibility condi- 
tions. There exist (5o,T > so that if 



(4.31) 



00 



then there exists a unique solution {u,p,r]) to (14. ip on [0,T] that satisfies the estimate 



(4.32) 



uit) 
rj{t) 



<VT+T 



00 



00 



for all t e [0,T]. 

Proof. The result without surface tension follows from Theorem 2.1 of [23]. In the case with 
surface tension, we proved a lower-regularity version of this result in Theorem 2.5 of [23]. The 
higher-regularity version we state here can be proved using a straightforward modification of our 
method in [23j, based on our new a priori energy estimates stated in Theorem 13.41 □ 



4.3. Unified data analysis. In order to prove our nonlinear instability result, we want to use the 
linear growing mode solutions constructed in Theorem 12.41 (cf. (|4.22|) ) to construct small initial 
data for the nonlinear problem, written in the unified equivalent perturbed operator form (I4.13p . 
Due to the nonstandard growth rate estimates stated in Theorem 12.71 (cf. (|4.21|) ) which force us to 
derive the nonlinear estimates in the higher-order regularity context (cf. (|4.23p ). we cannot simply 
set the initial data for the nonlinear problem to be a small constant times the linear growing modes. 
The reason for this is that the initial data for the nonlinear problem must satisfy certain nonlinear 
compatibility conditions in order for us to guarantee local existence in the space corresponding to 



norm 



II 00' 



which the linear growing mode solutions do not satisfy. 



To get around this obstacle, it is noticed that the nonlinear problem is slightly perturbed from 
the linearized problem and so their compatibility conditions for the small initial data should be close 
to each other. In the context of the nonlinear compressible Navier-Stokes-Poisson problem, Jang 
and Tice [12] used the implicit function theorem to produce a curve of small initial data satisfying 
the compatibility conditions for the nonlinear problem which are close to the linear growing modes. 
Such argument is general and is also suitable for our problem, so we may state this result for our 
setting without the proof. 

Proposition 4.5. Let u^,,r]^, be the growing mode solution stated in Proposition \4.S\ and assume 



the normalization 







(4.33) 



1 . Then there exists a number lq > and a family of initial data 
,2 fuit'T 



for I £ [0,io) so that the followings hold. 

1. UQ,r]Q satisfy the nonlinear compatibility conditions required by Theorem \4.4\ foi^ « solution to 
the nonlinear problem (j4.13p to exist in the norm ||H||oo- 

2. There exists a constant C > independent of l so that 
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and 
(4.35) 



00 



Proof. See the abstract argument before Lemma 5.3 of |12j . 



□ 

4.4. Proof of main theorems. With Propositions I4.2H4.31 Theorem 14.41 and Proposition 14.51 in 
hand, we can now present together the 

Proof of Theorems \1.1\ and \1.SX We shall prove the two main theorems together with the under- 
standing that when we write A we take A = A defined by (j2.32p for Theorem ll.il while for Theorem 
ll.2l we take A = A* defined by (j2.33p . In each case, we have that ^ < A < A as stated in Proposition 
14.21 It is this property that all we need in our unified proof. 

First, we restrict that Q < i < io < 9q, where lq is as small as in Proposition 14.51 and the value 
of ^0 will be fixed later to be sufficiently small. For < i < iO; we let tig and r/g be the initial data 
given in Proposition 14.51 By further restricting l we may use (j4.35p to verify that (14.311) holds. 



which then allows us to use Theorem 



(4.36) 



i=0 



to find 



solutions to the system (I4.13P with 



fl{L) 



By (|4.2ip and (|4.34p . we can estimate 



(4.37) 



u{i) 



< Ce 



At 



< Cie 



At 



for some constant Ci > independent of l. 

Let us now fix 5 > as small as in both Propositions 14.21 and 14.31 a-nd let C5 > be the constant 
appearing in Proposition 14.21 for this fixed choice of 6. We then define 6 = mm{6, nr^}- Denote 



(4.38) 



sup 

t 



tit) 



00 



<6)' and T* 



sup 

t 



< 2ie 



xt 



With io smah enough, (I433])-(|435]) and guarantee that T*, T** > 0. Then for all t < 

min{T'', T*, T**}, we deduce from the estimate (j4.23p of Proposition 14. 2^ the definitions of T* and 
T**, and that 

2 



<Cs 



00 



A 

+ 2 



(4.39) 



< 



< A 



Cs 



00 
u''{s) 

t 







U [S, 



ds 



00 







ds + Cs 

00 

u'{s) 



u''{s) 



ds + CsCl' 



00 



ds 



C5{2i) 
2A 



ds + C2i'^e^^K 



00 



for some constant C2 > independent of i. We may view ()4.39p as a differential inequality. Then 
Gronwall's lemma implies that 

X'Si^) ' < C2i''e^^' + C2i\^' f Xe^' ds 
Jl [t)/ 00 -^0 



(4.40) 
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We then deduce from Proposition 14.31 and (|4.40p that 



(4.41) 



















2 


< c 










00 



L e 



u''{s) 
rj''{s) 



+ C / e^(*-^) 
Jo 

^0 
2A — A 

Here we have used the fact that 2A — A > 0. 
Now we claim that 

(4.42) T = min{r\r*,r**} 

by fixing small enough, namely, setting 

/ S 1 1 



ds 



00 



(4.43) 



iQ = mm 



\2V2^'2(Ci+C3) 
Indeed, if T* = min{TS T*, T**}, then by (fOO]) . we have 



(4.44) 



u'{T*) 



00 



which is a contradiction to the definition of T* . If T** = min{TS T*, T**}, then by (jMID, (jOT]) . 
and the fact that A/2<A<A, we have that 



< ie 



AT" 



which is a contradiction to the definition of T**. Hence ()4.42p must hold, proving the claim. 
Now we again use ()4.4ip and ()4.37p to find that 





2 Ay. 



u{l) 



2 2AT" 



(4.45) 



u'{T') 



> 



Le 



XT'- / u 



u{i) 



(4.46) 



> ,eAT' _ ^2^^gAT' _ ^3,2g2AT' 

> ,e^T' _ ,2^^g2AT' _ C3,2e2AT' 

> ^0 - C'l^o ~ C'36'0 > 

This completes the proof of Theorem 11.11 and Theorem 11.21 in the unified setting. 



□ 
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Appendix A. Analytic tools 

A.l. Poisson extension. We will now define the appropriate Poisson integrals that allow us to 
extend r]±, defined on the surfaces Ti±, to functions defined on Q, with "good" boundedness. 

Suppose that S+ = x {1}, where := (27rLiT) x {2-KL2T). We define the Poisson integral 
in T2 X (-00, 1) by 

?e(L-iz)x(L2-iz) 
where for ^ G (L^^^) x {L:^'^Z) we have written 

(A.2) fiO = I /(x')-^^Ji=dx'. 



t2 2'K^JLiL2 

Here "— " stands for extending downward and "1" stands for extending at X3 = 1, etc. It is 
well-known that : ^ ij^+i/2(-p2 ^ (_oo, 1)) is a bounded linear operator for s > 0. 

However, if restricted to the domain fi, we can have the following improvements. 

Lemma A.l. Let V-^if be the Poisson integral of a function f that is either in H'^{T,j^) or 
iir'?^^/2(S_|_) for (7 G N = {0, 1, 2, . . . }, where we have written H'^{I]+) for the homogeneous Soholev 
space of order s. Then 

(A.3) l|V«P_,i/||o < ||/||^,-V2(T2) IIVP-.i/llo < ll/lli(T2)- 

Proof. See Lemma A.3 of [lOj. □ 

We extend 77+ to be defined on Q by 

(A.4) ??+(x',X3) =V+r]+{x',X3) := V-^iri+{x' ,x-i), for X3 < 1. 

Then Lemma lA. 1 1 implies in particular that if ry_|_ G for s > 0, then 77+ G H^(yi). 

Similarly, for S_ = x {0} we define the Poisson integral in x (— oo,0) by 

It is clear that V-^ has the same regularity properties as "P-,!- This allows us to extend rj- to be 
defined on However, we do not extend r/_ to the upper domain by the refiection since this 
will result in the discontinuity of the partial derivatives in 3:3 of the extension. For our purposes, we 
instead to do the extension through the following. Let < Aq < Ai < • • • < \m < 00 for m G N and 
define the (m+1) x (m + 1) Vandermonde matrix V{\q, Ai, . . . , Am) by V^(Ao, Ai, . . . , \m)ij = 
for i, j = 0, . . . ,m. It is well-known that the Vandermonde matrices are invertible, so we are free 
to let a = (qq, ai, . . . , a^)^ be the solution to 

(A.6) V{\o,Xi,. . . ,Xm)a = Qm, 

Qm = (1) 1, • • • 1 1)"^- Now we define the specialized Poisson integral in T^ x (0, 00) by 

It is easy to check that, due to (|X6]l . dlV+flf{x', 0) = dlV-flf{x', 0) for ah < / < m and hence 
(A.8) 9°P+,o/(a;',0) = 5"P_,o/(x', 0) Va G with < |a| < m. 
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These facts allow us to extend r/_ to be defined on by 

/ A fi\ - / / \ -T-, / t \ f 7^+ o^-(a;', 2:3), X3 > 

(A.9) r]_{x',x^)=V-ri^{x',x^):= { ' ) ^ ' 

i P_,o??-(2; ,2;3), X3 < 0. 

It is clear now that if G for < s < m, then f/_ G H^{Vl). Since we will only work 

with s lying in a finite interval, we may assume that m is sufficiently large in ()A.6P for G H^{il.) 
for all s in the interval. 

A. 2. Some inequalities. We will need some estimates of the product of functions in Sobolev 
spaces. 

Lemma A. 2. Let U denote a domain either of the form Q± or of the form T,±. 

(1) LetO<r<si<S2 be such that si > n/2. Let f e H'^{U), g G H'^{U). Then fg G H''{U) 
and 

(A.IO) < 11/11^., Mhs, . 

(2) Let < r < si < S2 be such that S2 > r + n/2. Let f G H^^iU), g G H^-^iU). Then 
fg G H''{U) and 

(A.ll) 11/511^,. < 11/11^,.! Mhs. ■ 

(3) Let < r < si < S2 be such that S2 > r + n/2. Let f G H-''{T,), g G H'^iJ:). Then 
fg G and 

(A.12) <II/ILJ|5IL,- 

Proof. These results are standard and may be derived, for example, by use of the Fourier charac- 
terization of the spaces. □ 

The estimates in (2) and (3) in Lemma IA.2I above are not the best possible estimates of that 
form. We now record one specific improvement that we will use for products in qH^{Q,). 

Lemma A. 3. Suppose that f G H'^{ft) with s > 3/2 and g G oH^{^l). Then fg G oH^{^}) with 
II/5II1 < ll/llslblli- Moreover, if g G {oH^n))* then fg G {oH^^))* via {fg,ip), := (gjif),. In 
this case we have the estimate \\fg\\[QH^)* ^ 11/11811511(0^1)* • 

Proof. See Lemma A. 3 of [23j. □ 
We will also use the following lemma. 

Lemma A. 4. The following hold. 

(1) Let f G g G H^Q), then fg G {oH\n))* and 

(A.13) WfgWi.mm' ^ ll/lloll5lli. 

(2) Let f G g G H^/'^{T.), then fg G H-^''^{T.) and 
(A.14) ll/5ll-i/2<ll/l|o||5l|o. 

Proof. See Lemma A. 4 of [23]. □ 

In the following lemma we let U denote a periodic domain of the form Q.± with flat upper 
boundary and lower boundary F;, which may not be flat. We now record some Poincare-type 
inequalities for such domains. 

Lemma A. 5. The following hold. 



(2) ll/llL2(r„) < \m\\L^i^u) for f G H\U) so that f = on F;. 

(3) ll/llo < ll/lli < llV/llo for all f G H\U) so that / = on F;. 
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Proof. See Appendix A. 4 of [10]. □ 

We will need the following version of Korn's inequality. 
Lemma A. 6. It holds that \\u\\i < \\Ou\\o for all u G oH^in). 

Proof. See Lemma 2.7 of [1]. □ 

A. 3. Elliptic estimates. Considering the two-phase stationary Stokes problem 

-fj.Au + Vp = inft 

div u = F"^ in $7 

(A.15) I - /i+D(n+))e3 = on S+ 

M = 0, l{pl - //B(u))e3l = -F- on S_ 

= on Sb, 

we have the following elliptic regularity theory. 

Lemma A.7. Let r > 2. // G H'-^{G), F^ G H'-^iG), F^ £ H'-^/'^{V), then the problem 
(jA.lSP admits a unique strong solution {u,p) G {qH^^G) riH^{G)) x H'''^^{G). Moreover, 

(A.16) + < \\F\^2 + \\F\-1 + \\F^\\r-Z/2- 

Proof. See |23l Theorem 3.1]. The proof follows by using the flatness of the interface S_ and 
applying Lemma lA.81 below. □ 

We let G denote a horizontal periodic slab and write F for its boundary (not necessarily flat), 
consisting of two smooth pieces Pi, We shall recall the classical regularity theory for the Stokes 
problem with Dirichlet boundary conditions on both Pi and P2. 

—p,Au + Vp = / in G 

div u = h in G 

u = ifi on Pi 

ti = on P2. 

The following records the regularity theory for this problem. 

Lemma A.8. Let r > 2. Let f £ H''-^{G), h G H'-^{G), G ^-^/^{Ti), 992 G H'-^/'^{T2) he 
given such that 

(A.18) I h= VI- v+ m-^, 

Jg Jfi Jt2 

then there exists unique u G H^'{G), p G H^^^{G) (up to constants) solving (jATTj)- Moreover, 

(A.19) \\u\\Hr{G) + I|Vp||h'-2(g) ^ ll/lk'-2{G) + \\h\\m-^(G) + Ikilln— i/2(ri) + Il<^2||//r-i/2(r2)- 

Proof See [151 El]. □ 

A. 4. Pressure estimate. By the following lemma, the pressure can be viewed as a Lagrange 
multiplier, which in turn gives an L^ estimate for the pressure. 

Lemma A. 9. // A G {oH^{^))* is such that A{v) = for all v G oH^{^}), then there exists a 
unique p G L^(il) so that 

(A.20) {p, div v) = A{v) for all v G qH^ {^) 

and we have the estimate 

(A.21) \\ph<\\M\iomm*- 

Proof See [H]. □ 



(A.17) 
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1 /2 

A. 5. Commutator estimate. Let J = (1 — A) ' with A the Laplace operator on T", n > 1. 
We define the commutator 

(A.22) {J'.n9 = J\f9)-fJ'9- 

We recall the following commutator estimate: 

Lemma A. 10. 

(A.23) \\\J\M\l^ < W'^fh^W-'gh^ + II^VIIl^II^IIl-. 

Proof. See \13\ Lemma XI] for the case R". The case of T" can be handled similarly. □ 
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